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Elastic systems driven in a disordered medium exhibit a depinning transition at zero temperature 
and a creep regime at finite temperature and slow drive /. We derive functional renormalization 
group equations which allow to describe in details the properties of the slowly moving states in both 
cases. Since they hold at finite velocity v, they allow to remedy some shortcomings of the previous 
approaches to zero temperature depinning. In particular, they enable us to derive the depinning law 
directly from the equation of motion, with no artificial prescription or additional physical assump- 
tions, such as a scaling relation among the exponents. Our approach provides a controlled framework 
to establish under which conditions the depinning regime is universal. It explicitly demonstrates 
that the random potential seen by a moving extended system evolves at large scale to a random field 
and yields a self-contained picture for the size of the avalanches associated with the deterministic 
motion. At finite temperature T > we find that the effective barriers grow with lenghtscale as the 
energy differences between neighboring metastable states, and demonstrate the resulting activated 
creep law v ~ exp (— C/~ M /Tj where the exponent n is obtained in a e = 4 — D expansion (D is 
the internal dimension of the interface). Our approach also provides quantitatively a new scenario 
for creep motion as it allows to identify several intermediate lengthscales. In particular, we unveil 
a novel "depinning-like" regime at scales larger than the activation scale, with avalanches spread- 
ing from the thermal nucleus scale up to the much larger correlation length Rv- We predict that 
Rv ~ T~ a /~ A diverges at small drive and temperature with exponents cr, A that we determine. 



I. INTRODUCTION 



Understanding the statics and dynamics of elastic sys- 
tems in a random environment is a longstanding prob- 
lem with important applications for a host of experi- 
mental systems. Such problems can be split into two 
broad categories: (tj propagating interfaces such as r&ag- 
netic domain wallsu, fluid invasion in porous mediae! or 
epitaxial growtho; (ii) pcriodio-systems such as vortex 
latticestLcharge density wavesu, or Wigner crystals of 
electrons^. In all these systems the basic physical in- 
gredients are identical: the elastic forces tend to keep 
the structure ordered (flat for an interface and periodic 
for lattices), whereas the impurities locally promote the 
wandering. From the competition between disorder and 
elasticity emerges a complicated energy landscape with 
many metastable states. This results in glassy properties 
such as hysteresis and history dependence of the static 
configuration. In the dynamics, one expects of course 
this competition to have important consequences on the 
response of the system to an externally applied force. 

To study the statics, the standard tools of statistical 
mechanics could be applied, leading to a good under- 
standing of the physical properties. Scaling arguments 
and simplified models showed that even in the limit of 
weak disorder, the equilibrium large scale properties of 
disordered elastic systems are governed by the presence 
of impurities. In particular, below four (internal) di- 
mensions, displacements grow unboundedly^ with the 
distance, resulting in rough interfaces and, loss of strict 
translational order in periodic structuresa. To go be- 
yond simple scaling arguments and obtain a more de- 



tailed description of the system is rather difficult and 
at present only main two methods, each with its own 
shortcomings, have been developped. The first one is 
to perform a perturbative renormalization group calcu- 
lation on the disorder, and is valid in 4 — e dimensions 
to first order in e. InJ-his functional renormalization 
group (FRG) approachffl, the whole correlation func- 
tion of the disorder is renormalized. The occurence of 
glassiness is signalled by a non-analyticity appearing at 
a finite lengtscale during the flow, specifically a cusp in 
the force correlator. This yields non trivial predictions 
for the roughness exponents of interfaced!!. Another ap- 
proach relies on the replica method to study either the 
mean field limit (i.e large number of components) or 
to perform a gaussian variational approximation of the 
physical modpl. Using this variational approach both 
for manifoldsli3 and for periodic systems! 1 n! 12 !, correla- 
tion functions and thermodynamic properties could be 
obtained. It confirms the existence of glassy properties, 
with energy fluctuations growing as L where 9 is a posi- 
tive exponent. To obtain the glass phase in this method, 
one must break the replica symmetry. At a qualitative 
level, this is in good agreement with the physical intu- 
ition of such systems as being composed of many low 
lying metastable states separated by high barriers. As 
was clearly shown in the case of periodic manifolds, the 
correlation functions can be obtained by both the FRG 
and variational approach and are found to be in very 
reasonable ragreement, bridging the gap between the two 
methodsliIrl!3. Taken together, these two approaches. thus 
provide a very coherent picture for the staticsOtJ. In 
particular they allow to understand that although disor- 
der leads to glassy features in both the manifold and the 
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periodic systems, these two types of problems belong to 
quite different universality classes in other respects-such 
as the large distance behaviour of the correlationstJ. 

These properties have drastic consequences for the dy- 
namics of driven systems in the case, important in prac- 
tice, where an elastic description holds (i.e when plas- 
tic deformations can be neglected). Determining the re- 
sponse to an externally applied force is not only an inter- 
esting theoretical question, but also one of the most im- 
portant experimental issues. Indeed in most systems the 
velocity v versus force / characteristics is directly mea- 
surable and is simply linked to the transport properties 
(voltage-current for vortices, current-voltage for CDW 
and Wigner crystals, velocity-applied magnetic field for 
magnetic domain walls). In the presence of disorder it is 
natural to expect that, at zero temperature, the system 
remains pinned and only polarizes under the action of a 
small applied force, i.e. moves until it locks on a local 
minimum of the tilted energy landscape. At larger drive, 
the system follows the force / and acquires a non-zero 
asymptotic velocity v. In the simplest cases, the effect of 
disorder at large velocity is washed out and one recovers 
the viscous flow, as in the pure case. In the thermody- 
namic limit, it is believed that there exists a threshold 
force f c separating both states, and that a dynamical 
transition occurs at f c called depinning, where the veloc- 
ity is continuously switched on, like an order parameter 
of a second order transition in an equilibrium systcmll3, 
leading to a v—f characteristics such as the one shown in 
Figure R]. 
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FIG. 1. Typical force-velocity characteristics, exhibiting 
pinning at T = with a threshold force f c and creep at T > 0, 
f < fc- At large drive, the system flows as if there were no 
disorder. 

An estinaate of f c can be obtained via scaling 
arguments^ or with a critej-Ae-ip for the breakdown of 
the large velocity expansionOlll. Beyond f c , if one de- 
scribes the depinning as a conventional dynamical critical 
phenomenon, the important quantities to determine are 
of course the depinning exponent /3 giving the velocity 
V ~ (/ — /c)' 3 and the dynamical exponent z which re- 
lates space and time as t ~ r z . 

An even more challenging question, and experimen- 



tally at least as relevant, is the response at finite temper- 
ature T > 0. In the most naive description, the system 
can now overcome barriers via thermal activation, lead- 
ing to a thermally assisted flowE3 and a linear response 
at small force of the form v ~ ez4Z. T /, where A is some 
typical barrier. It was realizec£3~L3 that because of the 
glassy nature of the static system, the motion is actu- 
ally dominated by barriers which diverge as the drive / 
goes to zero, and thus the flow formula with finite barri- 
ers is incorrect. Well below the threshold critical force, 
the barriers are very high and thus the motion, usually 
called "creep" is extremely slow. Scaling arguments, re- 
lying on strong assumptions such as the scaling of energy 
barriers and the use of statics properties to describe an 
out of equilibrium system, were used to infer the small / 
response. This led to a non linear response, characteris- 
tic of the creep regime, of the form v ~ exp(— C /~^/T) 
where fj, = (D — 2 + 2£ oq )/(2 — Ccq) and Ccq is the rough- 
ening exponent for the static Z?-dimensional system. 

Obtaining a detailed experimental confirmation of this 
behaviour is a non trivial feat, in reasons of the range 
in velocity required. Although in vortex systems these 
highly non linear flux creep behaviours have been mea- 
sured ubiquitously, it is rather difficult to obtain clean de- 
termination of the exponents, given the many regimes of 
lengthscales which characterize type II superconductors! 
In some recent measurement, some agreement with the 
creep law in the Bragg glass regime was obtainedo. 
Probably the most conclusive evidence for the above law 
was obtained, not in vortex systems, Uut for magnetic 
interfaces. Quite recently Lemerle et aln successfully fit- 
ted the force- velocity characteristics of a magnetic do- 
main wall driven on a random substrate by a stretched 
exponential form v ~ exp— /~ ' 25 over eleven decades in 
velocity. This provided evidence not only of the stretched 
exponential behavior, but of the validity of the exponent 
as well. 

Given the phenomelogical aspect of these predictions 
and the uncontrolled nature of the assumptions made, 
both for the creep and for the depinning, it is impor- 
tant to derive this behavior in a systematic way from the 
equation of motion. Less tools are available than for the 
statics, and averages over disorder should be made using 
dynamical methods. Fortunately, it is still possible to use 
a functional rcnormalization group (FRG) approach for 
the dynamical problem. Such-aanapproach has been used 
at T = to study depinningc3'E3. It allowed for a calcu- 
lation of the depinning exponents, in D = 4— e. However 
this approach is still rather unsatisfactory. The FRG flow 
used inE3'E£lis essentially the static one, the finite velocity 
being only invoked to remove - by hand - some ambigui- 
ties and to cutoff the flow, with no real controlled way to 
show that this is the correct procedure. Furthermore in 
these approaches it is also necessary to assume, instead 
of deriving them from the FRG, some scaling relations in 
order to obtain the exponents. Another rather problem- 
atic point is that, with no additional input, the method 
ofc3 would yield three universality classes for the depin- 
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ning: two universality classes depending on the nature 
of the disorder (random bond versus random field) for 
manifolds and one for periodic systems, while numerics 
and physical arguments^ suggested that only two (ran- 
dom field and periodic) universality classes could exists. 
In addition, since this is also intrinsically a T — (and 
v = 0) approach, it can not be used to tackle the creep 
behavior. 

We propose here a single theory for describing all the 
regimes of a moving elastic system, including depinning 
and the non-zero temperature regimes. Our FRG equa- 
tions contain from the start the finite velocity and finite 
temperature. They thus allow to address questions which 
are beyond the reach of either approximate scaling theo- 
ries, or v = FRG flow. For the depinning we are able 
to determine the conditions required for the existence of 
a universal depinning behavior, as well as computing the 
depinning exponents (and estimating / c ). We show in 
particular that only two universality classes exist (out of 
the three) for the depinning since we explicitly find that 
random bond systems flow to the random field univer- 
sality class. We can also extract from our equations the 
characteristic lengthscales of the depinning. The main 
advantage of our approach is of course to address the 
finite T small v regime as well. The method allows to 
derive the creep formula directly and thus allows to con- 
firm the assumptions made on the scaling of the energy 
barriers. In addition we show that the creep is followed 
by a depinning-like regime and determine its characteris- 
tic lengthscales. A short account of some of these results 
was presented in Rcf. |27|. 

The paper is organized as follows: in Section |n| we 
present the equation of motion and the types of disor- 
der studied here. Section III is devoted to a brief review 



of scaling arguments and a summary of useful results 
from perturbation theory, presented in Appendix |^. Sec- 
tion IV contains the field theoretical formulation of the 



problem and the associated renormalization group flow 
equations, derived in Appe ndix |^. The static case is 
studied in Subsection IV C , focusing on the appearance 
of the cusp. The effect of the temperature is studied in 
details in Appendices R and ||. In the next sections, we 
study the depinning dV) and creep (VI) regimes. Both 
sections contain the outline of the derivation and a phys- 
ical discussion. Appendix ^ is devoted to the effect of a 
small velocity on the FRG. We conclude in Section VII, 
referring to an extension of our work proposed in Ap- 
pendix In Appendix ^ we fix the notations used 
throughout the paper. 



II. ELASTICITY AND DISORDER 



Elastic systems are extended objects which "prefer" to 
be flat or well ordered. We are dealing with two different 
types of elastic systems which however can be treated in 
the same way. On the one hand, interfaces, i.e. surfaces 



with a stiffness that makes local distortions energetically 
expensive, on the other hand, lattices with elastic dis- 
placements allowed about a regularly ordered configura- 
tion. 

The first type is the easiest to visualize. The interface 
is assumed to have no overhangs and is thus described by 
a height function u r defined at each point r (see Figure^). 
Its energy is proportional to its area J yl~+ |Vu| 2 and 
in the elastic limit iVitl <C 1, reduces to 



H ol [u] 



|Vu| 



(2.1) 



relative to the flat u r — configuration (notations are de- 
fined in the Appendix |X|). We denote by c the stiffness, 
or elastic constant. 




FIG. 2. Top: an interface with height field u r above r. We 
denote by r the (internal) coordinates along the interface and 
by z the height coordinates. Bottom: a lattice with reference 
positions R and displacements ur from R. 

Periodic structures, such as flux line lattices or charge 
density waves (CDW), can be described by the same type 
of elastic Hamiltonian. For each point (or line) in the 
elastic periodic system one can introduce a (vector) dis- 
placement field ur that gives the shift from the reference 
position R (see Figure 0). The elastic energy for small 
displacements is given by a quadratic form in the dif- 
ferences ur — uri bet ween neighboring points and thus 
can be written as ( |2.l| ) in a continuum description (r be- 
ing a generic point in space). When u has more than 
one component, c should be understood as a tensor (see 
Appendix |F]) . 

To take the quenched disorder into account in such sys- 
tems it is necessary to express the energy of the above 
elastic structure in the presence of impurities. The cou- 
pling to a substrate or to local fields is easily written for 
interface models and is more subtle for lattices. Quite 
generally the coupling to disorder leads to an energy: 



Ha 



V(r,u r ) 



(2.2) 



which gives rise to a pinning force F(r,u) = —d u V(r,u) 
acting on the displacement u r . Depending on the micro- 
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scopic origin of the disorder term V, the coupling (|2 
leads to quite different phy sics . 

In the case of interfaces (2.2) originates from 



Ha 



p(r, z) 



V(r,z)p(r,z) 



iK z -(z — U r ) 



S(z — u r ) 



(2.3) 
(2.4) 



in terms of the density p(r,z). One then usually dis- 
tinguishes two cases: either "random bond" (RB) when 
V(r, z) is short range (random exchange for magnetic do- 
main walls), or "random field" (RF) as discussed below, 
where V(r, z) has long range correlations. 

In the case of periodic structures, the density p(r) can 
be expressed using the set of vectors k of the reciprocal 
lattice and (2.2) originates from 



H dis = / W(r)p(r) 

J r 



(2.5) 
(2.6) 



where po the average densityQ. The potential W is ran- 
dom, of short range Tf (e.g. point impurities for a vortex 
lattice or a CDW). We call this case "random periodic" 
(BP). 



In both cases, using (2.5 ,2.4, pT^j2.6| ) and (2.2) one ob- 
tains for the correlations of V in (|2.2[) 



(V(r,u) - V(r',u')) 2 = -2S rr >R(u - u') 



(2.7) 



where R(u) is a periodic function with—the periodicity a 
of the lattice in the periodic (RP) caseo. The S function 
is cutoff at the microscopic scale r/. 
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FIG. 3. up: RB case, down: RF case. Right: correlator of 
the potential, left: correlator of the force. 

For an interface, R(u) has the shape shown in Figure [|. 
In that case, the width 77 of R(u) is typically given by 
the width of the interface or the size of impurities. The 
force resulting from such a random bond disorder has 
correlations!^!! 



as shown in Figure [3] where 

A(u) = -R"(u) 



(2.9) 



The signature of such a RB disorder for the interface is 
that J A = since R'(u) decreases to zero at infinity. 

Another type of disorder occurs in the case of interfaces 
separating two phases, like e.g. a domain wall in a disor- 
dered magnet. A random field couples differently to the 
two phases on the right and left of the interface, thus the 
energy resulting from the coupling to disorder involves an 
integral in the bulk of the system and not just at the in- 
terface positio n. T he correlation of the force can still be 
expressed by (2.8) and A still decreases to zero above a 
scale Tf as shown on Figure 0. Contrarily to the RB case, 
J A does not vanish. For a single com ponent displace- 
ment field u, the RF, of correlator ( |2.8| ), is still formally 
the derivative of a potential V(r,u) = — J u du' F(r,u'). 
The correlations of this fictious potential are of the form 

(|3) with R(u) = - Jq'W Jo' du"A(u"), and one has 
R(u) ~ ~ J A for \u\ 3> Tf which can be visualized 
as a random walk (where u plays the role of "time" and 
the random field strength J A — — 2R'(oo) is the "diffu- 
sion constant"). Contrarily to the RB for which R(u) is 
short range, R(u) for the RF grows at large u as shown 
on Figure |3J. 

In this paper we study the overdamped driven motion 
of such elastic systems which obey 



rjdtUrt = c\7 2 u rt + F(r, u rt ) + ( rt + f 



(2.10) 



where r\ is a friction, / is the external driving force den- 
sity and £rt a Langevin noise. The correlation (CrtCr'f) — 
2r)T5 rr >Stt' defines as usual a temperature T for th is ou t 
of equilibrium system. The long time behavior of (2.1C) 
at zero drive / = 0, r educe s to the thermodynamics 
at temperature T. In ( 2.1 0| ) the bareE3 pinning force 
F{r , u) is gaussian with zero average and correlator given 
by ( |2.8| ). We will consider three universality classes for 
A corresponding to an interface in a random potential 
(RB), in a random field (RF) or a periodic system in a 
random potential (RP). Physical realizations of such dis- 
orders would be respectively a random anisotropy forj-a. 
magnetic domain walltl, the ptandom field Ising systemsEJ 
and vortex lattices or CDWlBr3. 

It is also useful to rewrite ( |2.10| ) in the co-moving frame 
at average velocity v — (dtu r t). In the remainder of this 
paper, we switch to u rt — * u rt + vt and thus study the 
following equation of motion 



{rjd t 



(dtUrt) 

cV 2 )u rt 





F(r, Vt + U rt ) +Crt+ f 



(2.11) 



F(r, u)F{r>, u 1 ) = 5„.A{u - v!) 



(2.8) 



where / = / — nv is the average pinning force and r be- 
longs to a D-dimensional internal space. From now on 
we specialize to an unidimensional displacement field u r 
as would be the case for an interface model or a single 
Q CDW. This simpler case already captures the main 
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physics at small velocity, investigated here. Extensions 
to many-component systems will be briefly discussed. 

Before giving a quantitative treatment using renormal- 
ization group, let us review the qualitative arguments 
which have been given previously to describe the physics 
originating from (2.10). 



III. PRELIMINARY ARGUMENTS 



1. Statics 



In the absence of drive, ( [2.10 ) is equivalent to the equi- 
librium problem at temperature T. The state of the sys- 
tem results from the competition between elasticity, pin- 
ning and thermal fluctuations. The physics of suet 
lems can be investigated by a host of methodaJEIEf 
and here we only recall the salient points. Temperature 
does not play an important role as will become clear and 
we begin with the T = case. 

A subsystem of size i?, with displacement w(R) = 



{ur — uq) , is submitted to a typical elastic force den- 
sity /ei = cw(R)/R 2 and to a typical pinning force den- 
sity /pi n = y/ A(0) JRP. Balancing these quantities, one 
obtains that elasticity wins at large scales for D > 4, 
resulting in a flat interface with a priori bounded dis- 
placements. In D < 4, systems of size R smaller than 
the Larkin length 



R,= 



2 2 \ D 



A(0) 



(3.1) 



wander as predicted by the Larkin mode 



10: 



w(R) ~ r f 



(3.2) 



At larger scales R > R c , the system wanders further 
than the correlation length rf of the disorder. This sim- 
ple picture breaks down and the system can be viewed as 
made of Larkin domains of size i? c , which are indepen- 
dently pinned. First order perturbation theory confirms 
this picture below the Larkin length. The static equilib- 
rium (equal-time) correlation function at T = is (see 
Appendix |b|) 



_ A(0) 

The wandering computed from (|3.3| 



-(it r ,t - Uo,t) 



(3.3) 



(3.4) 



for D = 4 — 6 gi ves back (3.2), and we recover the scaling 

2. 



Rc, 




(3.5) 



We used that / kzS^LL = A D r^ D for 2 < D < 4 with 

A D=4 - e = -W 2 - 2 r[-e/2]/16 - Si/e when e -► 0+ 

The remarkable feature is that T — perturbation the- 
ory (either using replicas or equilibrium dynamics) gives 
that (3.3) is exact to all order£3v3 in A,_and is identi- 
cal to the correlation in the Larkin model] . Indeed, the 
naive perturbation series organizes as if the pinning en- 
ergy were simply expanded in u (thus the pinning force 
is independent of u with F(r)F(r') = A(0)<5(r — r')), 
resulting in a gaussian model. 

In fact, due to the occurence of multiple minima 
beyond i? c , this perturbative result is incorrecttj at 
large scale. It can be shown, for example on discrete 
systems, that if a configuration u^ s which minimizes 
H[u] — f [§(Vu r ) 2 +V(r,u r )] is defined on a volume 

larger than R c , then the Hessian J M [u GS ] becomes 

smg ularE Such instability appears clearly in a func- 
tional ranormalization group (FRG) treatment of the 
problems which proves that A becomes nonanalytic be- 
yond the length R c , as will be discussed below. It can 
also be seen within variational or mean-field treatments 
using replicasliy that replica symmetry breaking (RSB) 
is necessary to describe the physics beyond the Larkin 
length R c . Using either replicas with RSB or the FRG 
it is possible to describe the physics at all scales and to 
obtain the correct roughness exponent £ eq defined by 



R 

w(R) ~ rf [ — 



(3.6) 



wheije the value of £ eq depends on the statics universality 
classo. Since disorder induces unbounded displacements, 
the system is rough and the temperature is always for- 
mally irrelevant in D > 2. It is described by a T = 
fixed point, characteristic of a glass phase. 



2. Depinning 

An elastic system does not necessarily move under 
the action of a driving force. The disorder leads to 
the existence of a threshold force f c at T = as 
shown on Figure FjL- .A simple dimensional estimate of 
f c can be obtainedE3 by computing the sum of the in- 
dependent pinning forces acting on the Larkin domains 
(R/ R C ) D y/AffijR^ and balancing it with the driving 
force acting on the same volume R D f. This gives 



f 211 
/c ~ Rl 



(3.7) 



expression (3.1) by equating the wandering to r?' 



Another estimate of f c comes from th e larg e velocity 
expansionEaO of the equation of motion ( 2.11 ) (from the 
criterion / ~ rjv). It coincides with (3.7). 
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F° r / ^ fc the system moves with a small velocity, and 
it has been proposedtJ that depinning can be described 
in the framework of standard critical phenomena, with 
the velocity as an order parameter. This leads . to t he, 
assumption of two independent critical exponentsc3c3'E^I 
£ and z, defined through the correlation function in the 
comoving frame (in the stationary state for / — * fj) 



U 



o) 2 = r 2 ^C{t/r z ) 



(3.8) 



C(x) — ► est for x — ► and C{x) ~ x 2 ^' z for x — » oo. The 
dynamical roughening exponent £ close to the threshold 
a priori differs from its equilibrium value £ cq . Several re- 
lated exponents can be also introduced such as: (i) the 
depinning exponent (3; (ii) the correlation length expo- 
nent v describing the divergence of the length £ defined 
from the equal-time velocity-velocity correlation func- 
tion. They satisfy 



V ~ (/ - fcf 
S ~ (/ - /c)-" 



(3.9) 
(3.10) 



NumericallyEZI the motion of the system looks like a de- 
terministic succession of avalanches of size £ with charac- 
teristic time r ~ (/ — f c )~ zu ■ Fr prn t he argument u ~ wt 
and the statistical tilt symmetryc3Ea (see below), the ex- 
ponents P and ^ are usually determined from £, z by the 
scaling relations 



1 







2 - C (* - 



(3.11) 



To obtain these exponents analytically, one needs to 
perform an FRG analysis of the equation of motion. This 
will be discussed in more details in Section [v|. 



3. Creep 

At finite temperature T > 0, motion occurs at any 
drive. For low temperatures and very small drive / -C f c 
one expects the motion to be very slow, and thus, al- 
though it is a dynamical problem, a qualitative under- 
standing can be obtained by considering thermal acti- 
vation over barriers determined from statics arguments. 
An original estimatetj of such barriers led to linear, al- 
beit activated, response. However the effects linked to 
the glassy nature. oL the problem were understood at a 
qualitative levelcJc^l using scaling arguments. 

The argument proceed as follows: systems larger than 

R c have a (static) roughness w(R) ~ Tf (-ff^ an d 
hence the energy has typical fluctuations of order 



E(R) 



D-2+2Co. 



(3.12) 



with U c = cR®~ 2 r 2 the energy scale of a Larkin do- 
main. Assuming that the energy landscape is character- 
ized by a unique energy scale, and thus that the energy 



differences between neighbouring metastable states is the 
same as the energy barrier separating them as schemat- 
ically shown in Figure |], one obtains that the barriers 
height scale with an exponent D — 2 + 2£ eq . 




configuration 

FIG. 4. Energy landscape, with many metastable states in 
the valleys, differing by AE, and barriers U between them. 

Since the motion is very slow, it is usually argued that 
the effect of the drive is just to tilt the energy landscape, 
and the effective barrier becomes 



D-2+2Co. 



(3.13) 



The maximum of (3.13), obtained at R, 



•-opt 



Rc yf~) ' gives via Arrhenius law the largest time 

spent in the valley by the thermally activated system and 
thus yields the velocity 



exp 



Ik (L 



,= ^±^ (3.14) 

^ Seq 



known as the creep motion, characterized by the 
stretched exponential with exponent /i. Note that the 
effective barrier given by the above formula vanishes at a 



scale Rq ~ R c I 4- ) which d iverges as fast as Ropt , 

the typical size of a thermally activated excitation. 



R D-2+2C,_j. R D+£ 




FIG. 5. Effective barrier and motion by nucleation. 

This elegant scaling argument leading to the creep for- 
mula relies however on strong assumptions and does not 
yield any information on the detailed behavior, in par- 
ticular on what happens after the thermal jumps. The 
fact that static barriers and valleys scale with the same 
exponent is already a non-trivial hypothesis about the 
structure of the infinite-dimensional energy landscape. 



G 



Refined simulationsESTii] of a directed polymer D = 1 in 
N = 1 and N — 2 are consistent with the "equal scaling 
assumption" for this particular case, but a general proof 
is still lacking. The second, and more delicate hypothesis 
is the validity of the Arrhcnius description: (i) the system 
being out of equilibrium, it is not clear that dynamical 
barriers can be determined purely from the statics; (ii) 
one assumes that the motion is dominated by a typical 
barrier. These assumptions can turn incorrect for some 
specific problems. For example, in the case of a point 
moving in a one-dimensional randpsn potential, the v-f 
characteristics at low drive is nov-i of Arrhcnius type. 
Although this + 1 case is peculiar since the particle has 
no freedom to pass aside impurities (it is dominated at 
T = by the highest slope of the potential and at finite 
T by the rare highest barriers), one should also address 
the question of the distribution of barriers in higher di- 
mensions. 



IV. DYNAMICAL ACTION AND 
RENORMALIZATION 

A. Formalism and exact relations 



Let us now study the equation of motion (2.11) using 
a full FRG treatment. This will enable us to describe the 
physics at all lengthscales and in particular the depinning 
and creep regime. 

A natural framework for computing perturbation 
ibrium systems is the dynamical 



theory ir 
formalism 




tion (2.10 


) using 



mal noise and disorder can safely be done and yields the 
simple ( "unshifted" ) action 



S uns (u,u) 



iu rt (i]d t -cV )u rt -rjT / iu rt iu rt 

rt J rt 



-f / iu rt 
Jrt 

- [ iu rt iu r t'A(u rt — u rt ') 

Jrtt' 



(4.1) 



Disorder and thermal averages (A[it]) = (^4[u])s unB of 
any observable A[u] can be computed with the weight 
e _Suns . Furthermore, response functions to an external 



perturbation h rt added to the right hand side of ( 2.10| ) 
are simply given by correlations with the response field: 
(A[u]iu rt ) — j£—{A[u]). It can be checked that causality 
is satisfied: (A[{ut>}t'>t,u]iu rt ) vanishes. In the time- 
continuum, the reponse to a perturbation at time t of 
an observable depending on u t is ill-defined. We choose 
Ito convention for the equation of motion, which ensures 
that equal time response functions, and hence any dia- 
gram occurring in perturbation theory containing a loop 
of response functions, vanish. The continuum field theory 
necessarily breaks down at small scales and it becomes 



necessary to cut off the integrals over the modes at large 
q, using a large wave vector A. A full summary of the 
notations can be found in Appendix [A|. 

It proves more convenien t to work in the comoving 
frame (i.e. with equation ( 2.11 )). The corresponding 
action is 

S(u,u) = / iu rt (rjd t - cV 2 )u rt - rjT / iu rt iu rt 

Jrt Jrt 

-f [ iurt (4.2) 

J rt 

-\ [ iu rt iu rt > A(u rt - U rt ' +v(t-t')) 

Jrtt' 



where the field u satisfies {d t u r t) — (d t u rt }s = 0. This 



condition fixes / = / — rjv in (4.2). This quantity is the 
(macroscopic) pinning force, since it shifts the viscous 
law / = rjv by the amount of /. 

Several exact relations can be derived directly from 



(4.2). For any static field h r (vanishing at infinity) 
1. 



S(u+ -V~ 2 h,u) = S{u,u) 



Performing the change of variable u — > u + -jV 2 h gives 
Du Du u rt e" s( "' fi) = 

n 7-1' I i ^w-2i \ — S(u,u)+ f i&rthr 

Du Du (u r t H — V hr) e 'Jrt 



(4.3) 



Applying -^-\h=o yields the exact relation 



T^qt — — n 
cq z 



where we denote by TZ r t the exact response function. 
This symmetry, known as statistical tilt symmetry, en- 
sures that the elasticity is not corrected during the renor- 
malization. 

Another important relation can be derived from 



(4.4) 



-T? (dtUrt) S uns = / d t (u rt iUr't')s uns 
<¥ Jr't' 



This leads to the identity between the macroscopic mo- 
bility and the slope of the v-f characteristics at any drive 
and any temperature: 



— v(f) = lim -iojll q =o,u 



(4.5) 



This exact result can also checked explicitely in the case 
of a particle moving in a one-dimensional environments. 

To extract the physical properties from the action ( |4.2| ) 
it is necessary to build a perturbatiiie approach in the 
disorder. A particularly simple casdlj occurs when the 
velocity is very large. In that case the disorder operator 
in the action can be formally replaced by 
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i / iu rt iu rt 'A(v(t - t')) 

^ JrW 



(4.6) 



since one may neglect the u rt — u r t' compared to v(t — t ). 
This trick suppresses the non-linearity and the remain- 
ing action is quadratic. Furthermore, at large veloc- 
ity, A(v(t - t')) can be replaced by S(v(t - t')) J A = 
-S(t — t') J A and the disorder operator transforms into 
a temperature operator (because it becomes local in time 
t = t'). The resulting action is the dynamical action as- 
sociated to the Edwards- Wilkinson equations describing 
the motion of an elastic system in a purely thermal noise 



with (v r tv r 



(4.7) 



2r](T + T ev! )5rr'5w , Langevin noise El of 



additional temperature T cw — ^— ■ . 

Note that at T = the results at large v coincide with 
the pertubative expansion in powers of the disorder. The 
equal-time correlation function in the driven system with 
force / crosses over from the static i Larkin behavior 

at small scale to a thermal 4j behavior at larger scale 



A(0) 



for q 2 > EL 

rw 



for q 2 < 



(4.8) 



with the same T ow , generated at lengthscales r 3> 



B. Renormalization 

We renormalize the theory using Wilson's momentum- 
shell method. As the cutoff A/ = Ae~' is reduced, corre- 
sponding to a growing microscopic scale Ri = e l j A in real 
space, the parameters of the effective action for slow fields 
(whose modes q arc smaller than A;) are computed by in- 
tegration over the fast part of the fields (whose modes q 
lie between A; and A). This iterative integration gives 
rise to flow equations, better expressed in terms of the 
reduced quantities 



Ai(u) 



T, 



SpAf 

2 eC i)2 



D 



(cA 



S D A r 
cAfeW 
rjiv 



T, 



(4.9) 



fo = f~ Vov 

where Sn is the surface of the unit sphere in D dimen- 
sions divided by (2ir) D . The exponent £ is for the mo- 
ment arbitrary and will be fixed later so that the reduced 
parameters flow next to appropriate fixed points. In one 
case (RB) we will need a Z-dependent £, and it is under- 
stood that everywhere the rescaling factors e^' (appear- 
ing e.g. in fl4.9Q) should then be replaced by exp L dl' Cv ■ 



The reduced quantities A,T are homogeneous to u 2 and 
A to u. The parameter A/, which plays a crucial role 
below, can simply be expressed as the following ratio 



vt(R) \{R) 



Su(R) 



'7 



(4.10) 



of the distance (along u) travelled by the center of mass 
of the interface during t(R) and the roughness 5u(R) = 
r f (R/R c ) c . We have defined t(R) = rj(R)R 2 /c as the 
characteristic relaxation time in the model renormalized 
up to scale R. 

The details of the renormalization procedure can be 
found in Appendix |^. The flow equations read: 

dA(u) = (e - 2C)A(«) + (uA'{u) + TA"(u) (4.11) 

A"(u) (A((s' - s)A) - A(u + (s' - ,s )A ) 

-A'(u-s'A)A'(u + sA) 
+A'((s' + s)A) (a'(u - s'X) - A'(u + sA)" 



s>0. 



'>0 



91nA = 2-C - / e- s sA"(s\) 

Js>0 

ainf = e-2-2C+ f e- s sXA'"(s\) 

Js>Q 

df = e-^ l cAl f e- s A'(sX) 

Js>0 

where e = 4 — D and d denotes ^ . 

This complicated set ofpequations require a few com- 
ments: (i) as for the staticsB it is necessary to renormalize 
the whole function A, instead of just keeping few cou- 
plings as in standard field theory, (ii) the elasticity c is 
not renormalized dc = due to the statistical tilt sym- 
metry; (iii) our equations correctly show that no tem- 
perature can be generated at v = since the fluctuation 
dissipation theorem holds at equilibrium. 

Setting both T = and v = in ( fill 




the simplified set of equations used in Refs. 
ting only v = also yields equations found in Rcf. 
But compared to the previous FRG approaches of the 
dcpinning transition, our equations correctly take into 
account the effect of the velocity on the flow itself (in- 
stead of being tceated simply as a cutoff as in Ref. Eq) . 
Other attemptala to incorporate velocity and tempera- 
ture in the FRG equations did not obtain the first equa- 
tion giving the renormalization of the disorder at T > 
and v > 0. To be able to tackle the full dynamical prob- 
lem and study the depinning and the creep regime, one 
cannot avoid keeping track of the velocity and of the tem- 
perature in the flow, as will become clear later, since they 
yield non-trivial effects which are unreachable by simple 
scaling arguments. 

Our flow equations allow in principle to compute the 
whole v-f characteristics at low temperature. In the fol- 
lowing we analyse them in the three regimes correspond- 
ing to the statics (v = 0), to the depinning at zero tem- 



8 



perature (T = 0, / ~ f c ) and the creep regime (T > 0, 
/~0). 



C. Statics: the cusp 

At zero velocity, our approach is a dynamical formula- 
tion of the equilibrium problem. It thus allows to recover 
the known results about the statics, avoiding the use of 
replicas. The standard derivation of the statics using the 
FRG consists in writing a replicated Hamilonian for the 
elastic system pinned in a random potential with correla- 

7 -25 D (r- 



tor (V(r, u) - V(r', u')) = ~26 D (r - r ') R ( u n "')■ 
averaging over V the replicated action readsa 



After 



Silt] = 



(4.12) 



where a, b are the n replica indices. Performing an FRG 
analysis of (4.12) yields for the flow of R and T (remark- 
ably independent of n) : 

dR{u) = (e - 4C)R(it) + CwR'(u) + TR"{u) (4.13) 



+ ^R» 2 - R"(0)R» 



d\nT = e-2-2( 



Ri(ue cl ) and Ti 



r, - r -w s ° A ? 



cA? 



with R I ( u ) = e -4Ci^ II 

which are the same redefinitions as ( |4.9|) w ith the corre- 
lator A o f the force is related to R by fl2.9|) . It is e asy to 
see that (4.13) coincides with our equations (4.11) when 
v = which read 



dA(u) = (e - 2<C)A(u) + (uA'(<it) + 

+f A"(u) + A"(u) (A(0) - A(u) 
ainf = e-2-2C 



(4.14) 



A» 2 



Thus the two methods give the same results for the static 
and equal-time physical quantities. The additional in- 
formation conveyed by the flow of the fr ictio n r\ i n th e 



dynamical formalism is discussed later in V A and VB 



The temperature in the static system is an irrelevant 
operator, since it decreases exponentially fast with I. One 
thus commonly restricts to the T = version of the above 
equations. In that case, as is obvious from the closed 
equation 



dA"(0) = eA"(0) - 3A"(0) 2 



(4.15) 



the curvature A" (0) < (see Figure ||) of the correlator, 
for any initial condition, blows up at a finite length scale 
for D < 4 



I, 



3K(0)| 



(4.16) 



which corresponds to 



-/A 



3S D \A»(0)\ 



2 ? 

cr f 



iA 



S D A o (0) 



(4.17) 



when approximating |Aq(Q)| by Ao(0)/r 2 . One thus re- 
covers the Larkin length (|3.5[ ) . The blowup of the curva- 
ture of A corresponds to the generation of a cusp singu- 
larity: A becomes non-analytic at the origin and acquires 
for I > l c a non-zero A'(0 + ) < 0. However, the flow equa- 
tion for the running non-analytic correlator still makes 
sense. The non-analyticity just signals the occurence of 
metastable states. A well-defined fixed point function 
R*(it) exists for each of the RB, RF, RP cases when a 
suitable £ is chosen. 

In the RP case, £ = C oq = so as to conserve the 
period a, and the fixed point is given bytil 



A* (ax) 



ea 
~~6 



x(l — x) 



(4.18) 



for x G [0,1). 

In the RF case, £ = ( cq = e/3 so as to conserve the RF 
strength J A and the fixed point is given byQ 



= y - 1 - lny 



(4.19) 



A*(u)/A*(0), x = u^e/(3A*(0)) and 



where y 

A*(0) ~ Q.he 1 ^^ A ) 2/3 (see W§ 



In the RB case, it has been shownEl by numerical inte- 
gration of the fixed point equation that ( = £ oq ~ 0.2083 e 
yields a physical fixed point, for which no analytical ex- 
pression is available. 

Despite the irrelevance of the temperature, this opera- 
tor has important transient effects during the flow, even 
if we are left asymptotically with the T = cuspy fixed 
point. It can be shown (see Appendix ^) that the tem- 
perature hinders the flow from becoming singular at a 
finite scale. The running correlator evolves smoothly to- 
wards its cuspy fixed point and remains analytic, as was 
also noticed in Ref. ^?]. As shown in Appendix |5|, the 
rounding due to temperature occurs in a boundary layer 
of width proportional to T around the origin. This is 
confirmed by the ex iste nce of a well-defined expansion 
in T (see Appendix D3). This effect is missed by sim- 
ple perturbation theory that would naively suggest that 

the rounding occurs on a width proportional to "vT. In- 
deed the correlation function is proportional to T and 
smoothes A by A K — + A K e~ Tfi . Although not crucial 
for the statics this rounding has drastic consequences for 
the creep as analysed in Section VI. 

Let us return on the differences between the stati c and 
dynamical formalisms. Within the static approach ( 4.13| ) 
in the T — > limit, despite the occurence of the cusp at 
l c , the RG equation for R; (it) still makes sense after l c and 
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flows to a fixed point controlled by e = 4 — Q_ However 
the physical meaning of the cusp is delicately. On the 
other hand, the use of the dynamical formalism allows to 
put T — from the beginning but adds to the problem 
a time dimension and the corresponding parameter, the 
friction rj. In this dynamical version of the problem, the 
cusp has strong physical consequences which are more 
immediate: after l c , the cusp generates infinite correc- 
tions to the friction. This feature marks the onset of a 
non-zero threshold force at scales larger than the Larkin 
length and signals that an infinite time is needed to go 
from one metastable state to another. Metastability thus 
appears very clearly in the dynamical formulation of the 
statics problem. 

A simple physical picture of the cusp in the statics at 
T = was given in Ref. [l3|. The renormalized poten- 
tial V Ien (r,u) at scales R > R c develops "shocks" (i.e. 
discontinuities of the force — d u V vcn {r 1 it) of typical mag- 
nitude fdisc(R) at random positions). Let us now extend 
this description to draw the link with the critical force 
and to include thermal effects. 

The force correlator for small u — v! is dominated by 
the configurations with a shock present between it and 



1 U sc (R) 2 ^\u-u'\ (4.20) 



(F ren (r, u) - F len (r, u')) 



where dp/du denotes the probability to find a shock 
between it and u + du. Identifying the rhs with 
R~ D \A[. cn (0 + )(u — u')\ one finds, using the rescalings 



(4.9), that the discontinuity in the force has the following 
scale dependence 



/disc(-R) 



' { .Rc 



"(2-C) 



EE ff(R) 



(4.21) 



and can thus be identified with an "effective critical 
force" f° s (R) at scale R, whi ch will play a role in the 
following (see Subsection |yiC|) . At R = R c , ff(R) re- 
duces to the true critical force f c . 

The renormalized problem at scale R being the one 
of an interface in a potential V ren (r,u) with the above 
characteristics, one can now easily understand the result 
that the cusp of A; (it) is rounded on a width Ti/x at 
T > 0. Extending the previous argument, one expects 
a rounding of a shock if the barrier between it and it' 
is of order T. Since near a shock the potential is linear 
of slope /disc(-R), the barrier is f&i ac {R)\u — u'\, and the 
thermal rounding should thus occur in a boundary layer 
of width it given by 



f disc (R)uR 



D 



(4.22) 



Using the rescalings (4.9), this is indeed equivalent to the 
expression Ti/x for the width of the boundary layer in 
rescaled variables found in Appendix [d[ 



V. DEPINNING 

At T = and v — ► 0, our flow equations give a self- 
contained picture of the depinning transition. Thanks 
to our formalism, the problem is reduced to the math- 
ematical study of (4.11), which although complicated, 
requires no additional physical assumptions. To focus on 
the depinning transition, we must analyze the solutions 
of these equations in the regime of small velocitiy where, 
using (l.£), A; = o is small. We will examine the various 
regimes in the RG flow keeping in mind that A; increases 
monotonically with I. 



The equations (4.11) involve averages over a range 
it ~ A; and thus one naturally expects that, at least at 
the beginning of the flow, A; (it) remains close to the 
v = solution. The two functions will differ in a bound- 
ary layer around it = of width denoted by pi . Although 
the precise form of the solution for |u| < pi (e.g. whether 
the cusp persists at v > 0) is very hard to obtain analyt- 
ically, fortunately most of our results will not depend on 
such details. As we discuss below, the main issue will be 
to decide whether p\ <C A; or not, which is a well-posed 
mathematical question. 

Let us start by analyzing the flow up to the Larkin 
scale l c of the statics, at which the cusp occurs and the 
corrections to the friction become singular in the v = 
flow. Here at v > one enters at l c a regime where A; 
is close to its fixed point (see Appendix |l|). Within the 
boundary layer, the effect of the velocity is to decrease 
the singularities of the statics. As shown in Appendix 
the blow-up of the curvature A"(0) is slowed down by 
the velocity as 

9A"(0) = eA"(0) - 3A"(0) 2 

-9A 2 A"(0)A iv (0) +C(A 4 ) 

and the same is true for the friction 



din?? = -A"(0) 



3A 2 A' 



(0) + o(\ 4 



If the blurring of the singularity results in a suppres- 
sion of the cusp, i.e. if A; remains analytic, one should 
wonder whether the v ~ flow can really remain close to 
A* since the convergence to the fixed point is crucially 
dependent on the existence of the non-analyticity and in 
particular on the term — A'(0 + ) 2 in the flow of A(0) in 
(4.14). A hint that A/ can stabilize for a while at v > 



is obtained by noting that one has (see ( |G2| )) 
<9A(0) = (e - 2<)A(0) 




which has indeed the correct sign to give the same effect. 

Hence it is natural to expect for I > l c that A; (it) 
has reached everywhere a fixed point form except in the 
boundary layer. The correction to the friction, crucial to 
determine the v-f characteristics, reads 
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di\n m =- e-'sA'/isXi 



(5.1) 



s>0 



and thus depends on the values of A; (it) for u ~ A;. 
To estimate this expression, one must know whether the 
width pi of the boundary layer is smaller than A; or not. 

To summarize these preliminary remarks, the flow in 
the Larkin regime I < l c is similar to the v = flow and 
A ; >; c is close to A* except for \u\ < pi. We will now an- 
alyze in details the flow for I > l c under the assumption 
that 



pi < A ( 



(5.2) 



As mentionned above, the validity of (5.2) can in prin- 
ciple be established by a mathematical or a n umerical 
analysis of our equations. It turns out that (5.2) leads to 
the most physically reasonable results. The alternative 
case will be discussed below. 



A. Derivation of the depinning law 



For I > / c , called the depinning regime, and relying on 



(5.2), the flow of rj becomes 

9 i ln ?7i ~-A*"(0 + ) 



(5.3) 



The friction is renormalized downwards with a non-trivial 
exponent — A*"(0+) = — ^ with £ = | for the RF case 
(see Ap pen dix E 1 ) and £ = for the RP case (see Ap- 
pendix [E 2| ). For the random bond one would naively 
take the static £ cq . However our flow equations show 
that during the Larkin regime, the form of the disorder 
correlator evolves to a RF, and thus C = § a l so m this 
case. This non trivial effect of the transformation for the 
dynamical properties of a RB into a RF is discussed in 
details in Section |VB| . 

Since A; keeps on growing in the depinning regime, 
the assumption that A; (it) can be replaced by A* will 
cease to be valid. This occurs when A; reaches the range 
rf(l) of A[(u), correlation length of the running disorder. 
This defines a scale ly = lnAi?y given by Xi v — r/(Zy). 
Above this scale, one enters a regime where the correc- 
tions due to disorder are simply was hed o ut by the ve- 
locity, since the integrals over s, s' in ( 4.11 ) average com- 
pletely over the details of A; (u) . One thus ent ers t he Ed- 
wards Wilkinson regime. Perturbation theory (4.8) shows 
that the interface is flat fpi these large scales for D > 2, 
the disorder leading onlyc3 to the effective temperature 
T 

± ew- 

The family of systems indexed by < I < oo have all 
the same velocity v and the same slope df/dv(v). How- 
ever they have lesser and lesser singular behavior f(v). 
We can thus iterate the FRG flow up to a point where 
the theory can solved perturbatively (e.g. above ly). For 
the depinning, one can simply use the fact that the renor- 
malized action at / = oo is gaussian and its friction -q^ is, 



from (4.5) equal to the slope df/dv of th e depinning char- 
acteristics. Using the flow of A; in (4.11), the expressions 



for \i v , A; c with (|4.9|) and (H2|) lead to 



A; v 



exp [(2 - 0f3(l v - Q] 



Vl c v 



with 



/3=l-A*"(0+)/(2-C) 



(5.4) 



(5.5) 



which will turn to be the depinning exponent and we have 
defined a characteristic force F c = crf/R 2 . Note that F r 



is not exactly the critical force f c . Solving (5.4) gives 



Ry 

T c 

Vly 

Via 



F c \ < 2 -o/3 

Vla V J 

Vi c v \ 3 
F,- 



(5.6) 
(5.7) 



Since the system at ly is nearly pure, one has r/i v m rjoc 
and, integrating over v the derivative — rj^ « r\\ v , 
one gets 



F c 



f-f(v = 0+) 







(5.8) 



which shows that the depinning is characterized by an 
exponent j3 and a pinning force f c = f(v = + ) (yet to 
be determined). 

The flow of // allows to fix the value of f c . Instead of 
just computing f c we also show that the integration of the 
flow off/ provides a second way to derive the depinning 
law ( |5.8| ). Indeed, as discussed below, in our formalism 
the term proportional to v which was problematic in the 
previous approachesEj cancels naturally. 

In the theory renormalized up to ly, the short scale 
cutoff is Ry and one can use fisrt order perturbation 



theory. One has (see (B7) and (Bll)) 



f lv = - A' lv (vt)R l 



(5.9) 



Since in the renormalize d th eory the disorder is close to 
A* (with the rescalings (4.9)), and the friction hidden in 
the response function is such that Xi v matches the range 
of A; v , the velocity disappears from (|5.9|) which gives 



(5.10) 



where A is some constant and the only w-dependent quan- 
tity is ly. To connect /; v to the initial parameters, 
one has to integrate the flow fi v — fo = f Q v dldifi. 



Expanding difi in (4.11) at small velocity and using 
A'(sA) = A'(0+) + sAA"(sA) + 0{\ 2 ) one recognizes in 
the second term the correction to rj. Thus for l c < I < ly 
one has 
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dji ~ e-^ l cA 2 A*'(0+) - vdi 



(5.11) 



where we dropped the sub-dominant terms in velocity. 
The integration of the flow gives 



flv ~ fo — ~fc I 1 



_ e -(2-C)(iv-Jc) 



li- 



no) 

(5.12) 

where we defined f c = cAge-^ 2 -^^- 1 A *' (0 + ) |/(2 — C) ■ In- 
jecting fo = f— rjov and (5.10), we note that quite re- 



markably the t]qv cancel each other. We are left with 
/ - f c ~ e -<-*-M lv - l '\Ae - f c ) + m v v (5.13) 



and <\5 



We already know (5.7) that rji v 

e lv ~ v < 2 -<)i3 thus both terms r.h.s of ( |5.13l ) scale like 
v 1 ^ . This leads to the following result to lowest order in 



2 _R2 ^ r 

RP 



— ca 

12 W 



1 - f RF 
1-| RP 
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(5.14) 

(5.15) 
(5.16) 



where we used the fact that £ = e/3 (RF or RB) or £ = 
(RP) and the link between v = |A*'(0+)| and r f (RF) or 
a (RP) stated in Appendix |Bj. In addition, we assumed 
that r)i c has a regular behavior when v — > 0, a non-trivial 



point which we discuss in Subsection V C 



B. Discussion 



The approach of the previous Subsection V A allows 
us to obtain the characteristics of the T = depinning. 
We extract the depinning exponent f3, the pinning force 
f c and the characteristic lengthscales from the equation 
of motion without any additional physical hypothesis or 
scaling relation. Although the depinning problem the ex- 
ponent (3 a pd t he critical force were determined in previ- 
ous studiesc3c3, our method is an improvement in several 
ways. 

To get the depinning exponent and critical force, two 
main derivations exist in the litterature. One of them 
extends the static FRG formalism to the oui-pf equilib- 
rium depinning problem at zero temperatures, using .an 
"expansion" around an unknown mean-field solutionEj. 
Instead of directly looking at the renormalizcd correlator 
of the disorder, the method obliges to deal with the time 
correlation of the force, C(v(t — t')) in Ref. ^5|. This pro- 
cedure does not allow for a precise enough calculation of 
the v-f characteristics to demons trate the c ancellation of 
of the t/qv term (in our equation (5.12,5.13)). In order to 
obtain a depinning exponent (3 different from its "mean- 
field" value j3 = 1, it is necessary in Ref. to neglect 



by hand in the small v limit a term proportionnal to v 
against a term proportionnal to v 1 ^ with f3 < 1. Our 
method, that directly uses averaging over the disorder 
and properly takes into account the velocity in the flow 
of the renormalizcd action, allows to show explicitly the 
needed cancellation. nr-i 

The other analytical studyE3'c3 of depinning does not 
consider the renormalization before the Larkin length 
and assumes that the singularity is fully developped be- 
yond this lengthscale. This amounts to take as a start- 
ing point the equation of motion at zero velocity with a 
cuspy correlator for the force, and the Larkin length as 
the microscopic cutoff. Since the anomalous exponent of 
the friction is — A*"(0) which is ill defined for a cuspy 
correlator, one is forced in this method to argue that it 
should be replaced by — A*"(0 + ) which is finite. This 
prescription and scaling relations linking the roughness, 
the depinning and the time exponent 2 — A*"(0 + ), al- 
lows to extract the depinning exponent. In our method, 
the ambiguities that existed in Ref. ^6] to write the flow 
of rj beyond l c when using the zero velocity equations, 
and the trick — > + becomes a well-defined mathe- 
matical property of our finite velocity RG equations: if 
721) is confirmed, our approach directly shows that the 
*"(0 + ) prescription is the correct one and allows to 
prove directly the scaling relations, instead of assuming 
them, to obtain the exponent. 

Furthermore the occurence of the asymptotic Edwards 
Wilkinson regime in Ref. ^6] has to be put by hand as 
a cut in the v — RG flow. The important correla- 
tion length Rv (denoted Ly in Ref. ^) at which this 
regime takes place is thus not well under control and has 
to be estimated from dimensional analysis. In our case 
the depinning regime is naturally cut when our A, which 
tells how fast the system runs on the disorder, reaches 
the range of the flowing correlator. The scale Ry at 
which it occurs, and above which the non-linearities are 
washed out, can clearly be identified with the correla- 
tion length of the moving interface (or more precisely, of 
the veloci ty- v elocity correlation). The physical interpre- 
tation of (pl\), i.e., 



R c 



(5.17) 



with z = 2 — ^£ is the following: Ry is the scale at which 
"avalanches" occur in the driven deterministic sytem. 
The motion proceeds in a succession of such processes, 
where pieces of interface of typical size Ry depin over a 
distance rf(Ry /R c ) during a time r f(Ry / R c ) /v. 

In addition to providing a clean derivation of the de- 
pinning exponents and of the critical force, our equations 
contain new physics that was unreachable by the previous 
methods. 

Although in principle one would expect three univer- 
sality classes (RF,RB,RP) for the depinaing exponent, it 
was conjectured by Narayan and FisherE3 that the rough- 
ness exponent of the system at the depinning transition 
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for RB or RF is equal to the roughness exponent of the 
static RF case, ( — e/3. This result cannot be obtained 
by the approach of Narayan and Fisher or that of Natter- 
mann et al. since these authors did not include the ve- 
locity in their RG analysis, and simply treated the small 
v limit as v = 0. On the contrary, our flow equations for 
the correlator shows directly that a RB disorder does in- 
deed evolve during the flow towards a RF disorder, leav- 
ing only two different universality classes (RF,RP) for 
the dynamics against three for the statics (RB,RF,RP). 
Such evolution is shown on Figure ^, where an initial RB 
becomes dynamically a RF. In Appendix ^| we show that 
the correction to J A, which measures the RF strength 
of the disorder, grows as 



d A = 2A 



A" 



0(X 4 ) 



where we have used £ = e/3. This ensures that a moving 
system, even at arbitrary small velocity, sees an effective 
random field at large scale. 



A(u) 



0.1 v 




FIG. 6. The evolution of a random bond to a random field 
correlator obtained by numerical integration of the flow. The 
initial condition I = of the flo w is a RB (A(u) shown as 
a full line on u > 0). Following (4.11), the running correla- 
tor transforms into a RF as shown on the snapshot of A;(m) 
near l c (dashed line), as can be seen by comparing with the 
characteristic shapes of RB/RF shown in Figure | 



a scaling function of u/Xi). However in that case, the 
exponents should differ from the standard ones (unless 
some hidden and rather my ster ious sum rule would fix 
the value of the integral in ( |5.l|) ). In the absence of an 
identified fixed point, it is not clear wh ethe r universality 
would hold. Again this crucial point (5_^2) can be defi- 
nitely answered by an appropriate integration of (4.11). 
Thus the present approach, which clearly takes v into ac- 
count, identifies a s fl5.2|) the condition under which the 
trick used in Refs. |25|J26| gives the correct exponents. 

Another intriguing point concerns the continuity be- 
tween the v = and the v — > problems. Indeed, to 



derive the depinning law ( 5.14 ) we have assumed that 
r]i c remains finite as v —* 0. However, we should recall 
that in the non-driven case (v = and / = 0) , rji diverges 
at l c and thus rji c — ocE3. If there is any continuity in 
the RG flow as v — > then r/i c — > oo in this limit. In 
that case the consequence would be (see ( |5.8[ )) a mod- 
ification of the exponent (3 — ► /3/(l — a) if rji c ~ v~ a 
(or weaker logarithmic multiplicative corrections). We 
would then find for the depinning a different result from 
the conventional one. Since we are unable to solve an- 
alytically accurately enough the equation for rj around 
l c , one should resort to a numerical solution of our flow 



equations (4.11) to resolve this question. Using (4.11) it 
is necessary to check that r\\ c does not diverge as v — > 
like a pow er o f v so as to recover the standard depinning 
exponent (5.5). The question is of particular importance 
since, if really a finite-scale behavior, occuring near i? c , 
would control the macroscopic asymptotic behavior, then 
again one could wonder whether universality would hold. 

Therefore, the description of depinning in terms of a 
standard critical phenomenon may be risky. Indeed as 
clearly appears in our FRG approach, since the fixed 
point at v = is characterized by a whole function A* 
(i.e. an infinite number of marginal directions in D = 4) 
rather than a single coupling constant (as in usual critical 
phenomena) the effect of an additional relevant pertur- 
bation, here the velocity, can be more complex due the 
feedback of v itself on the shape of the function during 
the flow. This is particularly clear in the RB case which 
dynamically tranforms into RF. 



C. Open questions 



VI. CREEP 



Our FRG equations prompt for several remarks and 
questions. In the previous sectio ns, we have examined 
the consequences of the property (5.2) and established in 
that case that the values of the exponents were the ones 
proposed in Refs. ^5|,^(|[36| Although we consider it as 
unlikely, we have not been able to rule out the possibility 
that either pi ~ A/ (or even worse, p\ > A;) and thus we 
should examine the consequences of a violation of prop- 
erty ( |5.2| ). If pi ~ A/, it is not excluded a priori that 
there exists another "fixed point" behavior (e.g. with 



We now deal with the non-zero temperature case. The 
system can jump over any energy barrier and overcome 
the pinning forces, thus it moves with v > for any 
drive / > and never gets pinned. Let us now show 
how our equations (4.11) allow to investigate the creep 
regime that occurs when the system moves very slowly 
with / <C f c , at low temperature. 
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A. Derivation from FRG 

As for the depinning, we are interested in infinitesimal 
velocities. The bare Ao is th us ve ry small. The main dif- 
ference compared to Section V A is that the temperature 
is now finite as well. The main effect of T is to round 
the cusp in the flow. Since we are interested in extremely 
small velocities, we will consider Ao as the smallest quan- 
tity to start with. A non-zero temperature gives thus rise 
to a new regime in the RG flow, where the rounding of 
the cusp is due to temperature and not to velocity. This 
leads to the following regimes in the FRG flow shown on 
Figure |?]. We will examine the various regimes in the RG 
flow keeping in mind that again, A; increases monotoni- 
cally with I. 

Ltirkin thermal sum ration <h'i>hming Edwards-Wilkinson 



X 



l T l d 



"c -s -I "d l V 

FIG. 7. Characteristic scales and regimes for creep motion. 

Just as in the previous case, we expect a Larkin regime 
for < I < l c with small corrections. Above l c the dis- 
order reaches a regime where scaling is imposed by the 
temperature. Indeed since A; c <JC Ti a jx one can forget 
about the velocity and the FRG equations are very sim- 
ilar to the v = and T > case. In Appendix ^ we 
show that the temperature rounds the cusp on a bound- 
ary la yer u ~ Tijx and we obtain the explicit scaling 
form (Dl) 



Aj(u) ~ A,(0) - fif{u X /fi) 
f(x) = Vl + x 2 - 1 
X=|A*'(0+)| 



(6.1) 



which in the statics holds at all scales larger than a scale 
of order l c . Here, because we focus on v — > 0, the scan- 
ning scale A; e is smaller than the width of the boundary 
layer, and the flow of the friction reads in this regime 



dln m ~ -Af(0) 



f, 



(6.2) 



The temperature being irrelevant by power counting, the 
initial flow of T is 

dlnf = -9 (6.3) 

since the anomalous correction to T vanish as A — > 0. 
Here and in the following, 8 = D — 2+2( cq denotes the en- 
ergy f luct uation exponent of the static problem. Together 
with (6.2) it shows that the friction grows extremely fast, 
like exp e 01 . This is the thermal regime where motion only 
occurs via thermal activation over barriers. The veloc- 
ity is so small that the center of mass motion is unim- 
portant and the temperature essentially flows as in the 
v — problem. We have determined the flow in its ini- 
tial stages, and we now determine the scale at which this 
behavior ceases to hold. 



The flow equation flOg ) for rji together with the scaling 
function (6.1) for A for u <~ Ti/x shows that ( |6.2| ) holds 
only until the new scale It = In ARt defined as 



A, 



Ti T /x 



(6.4) 



For I < It the temperature remains the main source of 
rounding of the cusp. Above that scale one must take 
the velocity into account. 

In fact, this simple picture is not complete since, be- 
fore reaching It another phenomenon occurs, leading to 
another lcngthscalc. In the thermal regime the correction 
to T due to diso rder competes with the simple exponen- 
tial decay and ([3J^) breaks down. This physically ex- 
presses that motion in a disordered landscape generates 
a thermal noise (provided some thermal noise is already 
present). Using |y]), one has dlnf w -6 + 6x 4 A 2 /!T 3 
at small A. Thus the correction to T reverts at a scale 
l s = \nAR s such that A s ~ ffj 2 /x 2 - Note that l s < l T . 
Above l s the temperature does not decrease any more 
due to heating by motion. One can show using ( 4.11 ) 
that T saturates and does not vary much until the scale 
It- We call this intermediate regime l s < I < It the 
saturation regime. We checked it using a numerical inte- 
gration of the flow in this regime with the scaling form 
of the disorder fl6.1| ) . Analytically, if we suppose that af- 
ter l s , the correction of T due to disorder dominates —6, 
then one would have in this regime an invariant of the 

flow di (rf — 6% 2 A^ w 0. If this were true, it is clear- 
that the flow could never realize the condition A/ ~ Ti/x, 
possibility that is excluded on physical basis and by the 
numerics shown in Figure @. 
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x=l/T 



-0.6 ~0.8 r— . 

FIG. 8. The integration of the flow using ( p.l\) after l c and 

reduced variables xi = x^t/Ti and yi = x*Af /Tf. The dotted 
lines indicate the set of points where dixi = or diyi = 0. 
Some trajectories are displayed, with an arrow showing the 
direction of growing I. The initial conditions for creep are 
close to the origin, and closer to the a;-axis. 



Despite the saturation of the temperature, (6.2) re- 
mains true after l s . Thus the friction and A keep on 
growing and one finally reaches the scale It at which 
the scanning length A; crosses the boundary layer width 
Ti/X- 
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Above It, a rigorous analytical analysis of (4.11) be- 
comes difficult. We however expect, since the velocity 
controls now the boundary layer, a regime similar to the 
depinning regime at T = to occur. Using the same 
arguments than for the depinning, one obtains in that 
regime 



Thus 



<91n77 = -A*"(0 + ) 
dlnf = 2- D -2( 



(6.5) 
(6.6) 



leading again to a decrease of the temperature, even 
slightly accelerated by a negative 0(e 4 / 3 ) exponent. Let 
us call Id the depinning scale at which one enters such a 
depinning regime. From the above discussion it is very 
reasonable to expect that one goes directly from the sat- 
uration to the depinning regime, i.e. Id — It ^ est. How- 
ever we cannot strictly rule out the possibility of an in- 
termediate regime (divergent Id — It when v — > 0) during 
which the correction to the friction goes smoothly from 
positive (as in the thermal and saturation regimes) to 
negative values (depinning regime). Again, it would be 
useful to settle this point through a numerical solution of 
our flow equations. Note that in the RF a nd R P cases, 
the exponent £ and the fixed point A* in ( |6.5| ) are the 
same as in the statics. However in the RB case, one have 
used a ^-dependent £ which crosses over between £ eq for 
I < It and ( = e/3 for I > l,i corresponding to the change 
from RB to RF fixed points A*. 

In the depinning regime, motion now proceeds in a 
similar way than for the one studied in Section Here 
again at large enough scale, velocity will wash out the 
disorder for I > ly with ly determined by A; ~ r /(0- 
One then enters the Edwards- Wilkinson regi me. 

Let us now compute from the flow ( 4.11] ) the 
lenghtscales defined above (see Figure 0) . 

In the thermal regime l c < I < l s one can compute 
either by integrating its flow or by equating the boundary 
values to their expression. This gives 



exp [(2 - 
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where we defined f c = e-^f and U c = e R^-^. Express- 
ing the scales as a function of the velocity leads to 
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(6.8) 



With 11 = 9/ (2- Coq)- 

In the saturation regime l s < I < It we proceed in the 
same manner and obtain 
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(6.9) 
(6.10) 



Assuming Id ~ It, the depinning regime Id < I < ly 
follows directly and 

\ W ^ f exp [(2 - Q{3{l v - Id)} 



lUjLpS(h-lc) 



leads to 
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(6.11) 

i 

(6.12) 

the depinning exponent (and C, the 



with (3 

dynamical roughness exponent). 

We are now in a position to compute the character- 
istics f(v). We fix a small velovity v and solve the flow 
equations for A; , A; and Ti up to ly . This allows to relate 
fi v to the unknown fo- We can now use the fact that at 
the scale ly, the disorder is essentially washed out and 
a perturbative calculation of fi v m — is possible. 
Solving backwards we determine fo , wich is simply / — rjv 
where / is the real force applied on the system and rj = rjo 
the bare friction. 

The correction to / can not be neglected during the 
depinning regime, thus, using fo = f — i]o v , foe = and 
expressing J Q dl difi one has 



dldifi 



ch lx (2- 

2-C 



(2-Ceq)'d 



(6.13) 



In the thermal regime th ere is essentially no correction 
to the flow of /. Thus ( 6.13 ) is controled by the de- 
pinning regime and one should integrate essentially be- 
tween Id and ly. In fact due to the e xpon entially de- 
creasing behavior of the integrand in ( 6.13| ) the whole 
integral depends in fact only of the behavior at the scale 
Id- Assuming that Id ~ It, using fl6.7|]6.9]), one sees that 



-(2-Ccq)'d <g- v f or 



rjv 

— w exp 

Jc 



M = 



and thus one obtains 

(L 

T \f c 
D - 2 + 2Cc 



2-C 



(6.14) 



(6.15) 



eq 



The prefactor in front of the exponential cannot be ob- 
tained reliably at this order. Note that for the creep, 
contrarily to the depinning, the possible divergence of r)i c 
when v — > (and T — > 0) does not affect the argument 
of the exponential but only the prefactor. 
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B. Alternative method and open questions 

For the depinning it was possible to recover the de- 
pinning law using both the integration of t he flow of / 
and of the friction 77 and the relation (4.5). Although 
one can also use in principle this method for the creep 
it gives poor results in this case. Indeed contrarily to 
the derivation involving / one needs here the flow of 
77 in all regimes including the depinning regime I > Id, 
where r\ is still renormalized. Since the renormalization 
of 77 goes from large positive growth (first like expe el , 
then exponentially) in the thermal/saturation regime to 
negative in the depinning regime (where the system ac- 
celerates with subdiffusive z < 2) a precise knowledge 
of the behavior around It would be needed. Unfor- 
tunately the lack of precise analytical methods avail- 
able above Rt prevents from computing precisely such 
a crossover. A crude estimate of the flow can thus only 
give a bound of the exact result. If we use (e.g. in 
the RF or RP cases) the estimates of each regime, and 
the perturbative estimate of r\\ v in the theory at ly- 
Voo ~ Vlv + J t A'{ v {vt)tR& ~ e -(2-0av-ic) e ( it will ap _ 

pear that r\i v diverges faster than e~( 2 ~^' v ~' c - 1 when 



»-»0). The product of joTf j6.10| , |6l2| ) is equal to 



Integrated from to v, it yields 



fc 



exp 




(6.16) 



One would thus find, using the r\ method a non-Arrhcnius 
law for the creep regime. Even if one cannot strictly 
speaking exclude this result, as discussed above it is most 
likely an artefact of the approximate integration of the 
flow, and only a lower bound of the barrier height. In- 
deed compared to the integration of the flow of /, this 
procedure is much more sensitive to the neglect of the 
crossover It < I < Id- A more precise integration of the 
flow would very likely show a compensation between the 
latent growth of the friction during the decrease of d In rj 
(for It < I < Id) and the reduction of the friction occur- 
ing in the depinning regime Id < I < W ■ Note that if M- 
werc equal to r\i T then, one would recover (6.14). It would 
be useful to check explicitely on a numerical integration 
of the flow that such a cancellation does occ ur an d verify 
that the 77 method confirms also the result (6.14). 

We also note that the precise determination of the 
lcnghtscales for R > Rt depend on obtaining an accu- 
rate solution of the RG flow equations. In t he pre viou s 
section, we have obtained the formulas (6J) and ( 6.11 ) 
under some assumptions about the mathematical form 
of the solutions of the flow in the region where A; and 
Ti cross. These assumptions, discussed in the previous 
Section, should be checked further, e.g via numerical in- 
tegration. Although this should not affect the creep ex- 
ponent derived above, the precise determination of these 



length scales is important to ascertain the exact value of 
the scale Ry (i.e the avalanche scale discussed below). 



C. Discussion 



Since our flow equations ( 4.11 ) include finite temper- 
ature and velocity, they allow for the first time to treat 
the re gime of slow motion at finite temperature, directly 
from ( |2.10 ) . As for the depinning we derive directly from 
the equation of motion the force-velocity law and we ob- 
tain new physics. 

The first important result is of course the creep formula 
itself ( |6.14 ). Our method allows to prove t he ma in phys- 
ical assumptions, reviewed in Subsection III 3, needed 
for the phenomelogical estimate, namely: (i) the equal 
scaling of the barriers and the valleys; (ii) the fact that 
velocity is dominated by activation over the barriers cor- 
rectly described by an Arrhenius law. In our derivation 
such law comes directly from the integration of the flow 
equations in the thermal regime; (hi) the fact that one 
can use the static exponents in the calculation of the bar- 
riers. This appears directly in the formula (6.14) but can 
also be seen from the fact that in the thermal regime the 
velocity can essentially be ignored in the flow equations. 
We also recover the characteristic lengthscale predicted 
by the phenomenological estimate. Indeed one can iden- 
tify the scale (|j[ I6J) R T ~ i? c (/// c ) _1/(2_Ccq) as the 



i? op t of Subsection pTT 3 



Our equations allow to obtain additional physics in the 
very slow velocity regime. In particular, we see that the 
slow motion consists in two separate regimes. At small 
lengthscales R < Rt the motion is controlled by thermal 
activation over barriers as would occur at V = 0. This is 
the regime described by the phenomenological theory of 
the creep. Qualitatively, the main novel result obtained 
here is that the thermally activated regime is followed 
by a depinning regime, as shown by our equations. This 
leads to the following physical picture: at the length Rt, 
bundles can depin through thermal activation. When 
they depin they start an avalanche like process, reminis- 
cent of the T = depinning, up to a scale Ry. The 
propagation of the avalanche proceeds on larger scales in 
a deterministic way. Thus one is left with a depinning- 
like motion, and the size of the avalanches is determined 
by the natural cut of the RG (A = 77), i.e., the scale 
at which the propagating avalanche motion is overcome 
by the regular motion of the center of mass. One recov- 
ers qualitatively and quantitatively some features of the 
T = case at intermediate scale. The typical nucleus 
jumps over an energy barrier Uh ~ U c (Rt / R c ) e result- 
ing in v ~ exp — Hf- (^-jr^J . This jump of a region of size 

Rt initiates an avalanche spreadi ng over a m uch larger 
size Ry which we find to be (see (6.7,6.9,3.11)) 



Rv 
~R~c 



Rt 
~R C 



(6.17) 
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with v = (3/{z - C) = 1/(2 - C) and z = 2 - (e - C)/3 
the critical exponents of the depinning, and 8 the en- 
ergy exponent of the statics. Note that the correlation 
length Ry diverges at small drive and temperature as 
Ry ~ T-°f- X with a = % = ^ and A = ^- + ^. 

To push the analogy further one can consider that the 
avalanches at lengthscales R > Rt are similar to the 
ones occuring in a regular T = depinning phenomenon 
due to an excess driving force (/ — f c ) B g- Considering 
a minimal block size Rt instead of R c for this "creepy" 
depinning, Ry / Rt ~ (/ — /c) e fi > one obtains for this 
effective excess force: 



f; 



off 



/c 



T 



(6.18) 



linking the creepy motion at T > and the threshold 
depinning at T = 0. As explained before, there might 
be an uncertainty in the value of the avalanche expo- 
nent, wh ich c ould be changed by a quantity of 0(e). To 
confirm ( |6.17 ), one would need to further check the pre- 
cise behaviour of the solution of the RG equations for 
R > R T - 

One can understand qualitatively that the probl em a t 
scale R > Rt looks like depinning a ccord ing to ( |6.18 ). 
The tilted barrier (see Subsection [II 3| ) E{R, f) = 
U C (R/ R c ) e — f R D r f (R/ R c )^ c i to be overcome in order 
to move a region of size R (all barriers corresponding 
to smaller scales having been eliminated), vanishes atEll 
Ro(f) ~ Rt- For the T = depinning problem, one 
can define a scale dependent effective threshold force 
ff{R) ~ f c (R/R c )^ 2 -^ such that E(R,ff) = 
(also defined in Subsection [V C| ), which corresponds to 
the force needed to depin scales larger than R (the true 
threshold f c = f° s (R c ) being controlled in that ca se by 
the Larkin length). A possible scaling derivation of ( 6. IS ) 
is obtained by noting that at T > 0, non-activated mo- 
tion at scale R occurs when the tilted barrier E(R, /) is 
of the order of T. This yields a T dependent effective 
threshold force such that 



ff{R)-ff(R,T) 
ff{R) 



T 



U C {R/R c ) e 



(6.19) 



At R = Ro(f ), one has / = ff(R) and fl6T9|) is iden- 
tical to Q6.18 ) to zeroth order in e (i.e. (3 = 1). In fact, 
to apply the above static barrier argument, it might be 
better to work in the co-moving frame where the velocity 
of the interface vanishes. This amounts to replace / by 
/ — r]V in the previous argument, and E(R a (f), f) = 0, 
E(R {f)J — T}v) = T gives back ( fig ). 

The crossover between thermally activated processes 
and depinning-like motion can also be recovered by not- 
ing that the condition A; ~ Ti/\ which appears in the 
FRG flow can be rewritten as (using (fOj) and (4.10)): 



ff{R)vr{R)R L 



T 



(6.20) 



where the lhs is a natural energy scale involved in the de- 
pinning due to driving effect of the center of mass. If it 



is much larger than T, depinning effects dominate, while 
if it is smaller, the dynamics is activated. 

Finally many open questions still remain. Technically 
it would be interesting to reconciliate the two methods 
based on rj and / which proved to be equivalent for the 
study of depinning. In fact although the two methods 
should formally agree, the comparison at a given order 
in the RG is more subtle. Indeed -j^SJ = —6r], by integra- 
tion over I between and oo and derivation with respect 
tow, gives back rj^ = 4jj provided that /oo = 0. However, 

one should notice that in - d \5f — ~Srj, the derivative is 
understood at fixed parameters at the given scale. The 
occurence of this hidden dependence in the velocity in 
the running parameters makes the equivalence between 
both approaches delicate. However the additional term is 
of higher order in disorder. Thus, as pointed out above, 
it is very likely that a careful integration of the flow of 
r\ should resolve this discrepancy, but this remains to be 
explicitely checked. 

As for the T = depinning, the existence of the de- 
pinning regime at Id depends on the precise form of the 
boundary layer in the presence of a velocity. Note that 
the alternative scenario discussed in Section V_C, e.g. 
whether or not the depinnig regime is universal, would 
not affect the creep exponents, but only the subleading 
corrections. 



VII. CONCLUSION 

We examine in this paper the dynamics of disordered 
elastic systems such as interfaces or periodic structures, 
driven by an external force. We take into account both 
the effect of a finite temperature and of a finite velocity 
to derive the general renormalization group equations de- 
scribing such systems. We extract the main features of 
the analytical solution to these equations both in the case 
of the T = depinning (shown on Figure |l|) and in the 
"creep" regime (small applied force / and finite temper- 
ature). 

Our RG equations, when properly analyzed, allow to 
recover the depinning law v ~ {f—fc) 13 and the depinning 
exponent /3 also obtained by other methods. However, 
contrarily to previous approaches that needed additional 
physical assumptions, such as scaling relations among ex- 
ponents or by hand regularisation, our approach is self- 
contained, all quantities being derived directly from the 
equation of motion. It thus provides a coherent frame- 
work to solve the difficulties and a mbig uities encountered 
in the previous analytical studiesriri. In addition our 
method allows to establish the universality classes for 
driven systems. It shows explicitely that a random bond 
type disorder gives rise close to a random field critical 
behavior at the depinning. Thus the dynamics is char- 
acterized by only two universality classes (random field 
(RF) for interfaces and random periodic (RP) for peri- 
odic systems) instead of three. Since this phenomenon is 
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an intrinsically dynamical one, it was out of the reach of 
the previous analytical approaches that used v = flow 
equations together with additional physical prescriptions 
using e.g. the velocity as a cutoff on the v = RG flow. 

Of course one of the great advantages of the present 
set of RG equations is to allow for the precise study of 
the small applied force regime at finite T, for which up to 
now, only phenomenological scaling arguments could be 
given. Our FRG study confirms the existence of a creep 
law at small applied force 



rjv 

— exp 

Jc 



Uc 

T 



(7.1) 



with a creep exponent related to the static ones /z = 
(D — 2 + 2£ cq ) /(2 — Ceq)j where £ eq is the statics roughen- 
ing exponent. It provides a framework to demonstrate, 
directly from the equation of motion, the main assump- 
tions used in the phenomenological scaling derivation of 
the creep namely: (i) the existence of a single scaling for 
both the barriers and the minima of the energy landscape 
of the disordered system; (ii) the fact that the motion 
is characterized by an activation (Arrhenius) law over a 
typical barrier. 




FIG. 9. Schematic picture of the creep process emerging 
from the present study: while thermally activated motion oc- 
curs between scales R c (Larkin length) and Rt (thermal nu- 
cleus size), depinning-like motion occurs up to the avalanche 
size Rv- 

In addition, our study unveils a novel "depinning- 
like regime" within the creep phenomena, not addressed 
previously, even at the qualitative level since the phe- 
nomenological creep arguments did not address what 
happens after the thermally activated jump of the op- 
timal nucleus. Although the velocity is dominated by 
the time spent to thermally jump over the barriers, our 



equations show that the small / behavior consists in 
fact of two different regimes. Up to a size Rt motion 
can only occur through thermal activation over barri- 
ers. This is the regime described by the phenomenolog- 
ical approach to the creep. The optimal nucleus of the 
scaling estimate is given directly by the RG derivation 
as Rt ~ (l//') 1 /( 2 ~^ cq ). Remarkably, another interest- 
ing regime exists above this lenghtscale (see Figure ||). 
It emerges directly from our RG equations and can be 
given the following simple physical interpretation. In 
some regions of the system, bundles of size Rt depin 
due to thermal activation. These small events then trig- 
ger much larger ones, and the motion above Rt pro- 
ceeds in a deterministic way, much as the T — de- 
pinning. In particular once the initial bundle depins it 
triggers an avalanche up to a size Rv which is given by 
R v /Rt ~ (U c /T)»/P{R T /Rc) ei// P where 6, and v are 
the energy, depinning and correlation length exponents 
respectively. 

The present study also raises several interesting ques- 
tions which deserve further investigation, some of them 
rely on being able to obtain a more accurate solution of 
our flow equations. We have shown explicitly how to re- 
cover from our equations the conventional depinning law 
(and the scaling creep exponents). It rested on a math- 
ematical property, likely to hold, but not yet rigorously 
established, of the solution for the flow of the correla- 
tor of the disorder. Such behavior should be checked in 
details. The equations being quite complicated, a numer- 
ical solution, albeit delicate, seems to be appropriate. If 
the constraint (15.21) on the flow defined in Section |VB 



were found to be violated, then the conventional picture 
of the depinning would very likely fail, as we have ana- 
lyzed in detail. A similar question arises concerning the 
flow of the friction 77 as discussed in Section VIC. If 



the solution of the flow is found to depend on the pre- 
cise behavior at the Larkin length R c , it is likely that 
even universality could be questioned. These issues are 
a priori less important for the first, thermally activated, 
part of the creep regime, but because of the existence 
of a second, depinning-like regime, they would also have 
consequences for creep. Again, these question depend on 
the precise form of the flow and can be answered unam- 
biguously by a detailed enough analysis of our equations. 
It would also be of great interest to develop a more de- 
tailed physical picture of the crossover between thermally 
activated and depinning like motion since we found that 
both occur within the creep phenomenon. 

Several applications and extensions of our work can 
be envisioned. First, extensions to many-dimensional 
displacement field (of dimension N > 1), given in Ap- 
pendix [|], would be interesting to study within the meth- 
ods used here. One could check whether the approxi- 
mation used in Ref. ^ yields the correct result for the 
N > 1 depinning. Second, the effect of additional KPZ 
non-linearities could be investigated. In particular one 
could check the usual argument which yields that KPZ 
terms are unimportant for the depinningtH since their 



18 



coupling constant is proportionnal to the (small) veloc- 
ity. Also, extensions to other types of disorder, such as 
correlated disordeiH are possible. Finally, it should allow 
to describe in a systematic way the the thermal rounding 
of the dcpinning, i.e. the study of the v-f characteris- 
tics for / close to the threshold and small T. If one 
assumes that one can simply carry naive perturbation 
theory in T around the T = solution of the RG flow 
near f c (i.e. o nly k eeping the contribution beyond ly), 
one is led in (5.13) to an additional term proportional 
to T/v 2 , which readily yields the value for the thermal 
rounding exponent p = 1 + 2(3 proposed in Ref. [l8| (i.e. 
a scaling form near / = f c and small T for the velocity 
v T^/p$ ^Zz^k^). Although this exponent seems to be 

consistent with starting values A -C T, its validity could 
be further checked by solving our RG flow equations at 
small T. 
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APPENDIX A: NOTATIONS 



Here are some notations and conventions and diagram- 
matics we use in the text. 

The surface of the unit sphere in D dimensions di- 
vided by (2tt) d is denoted by S D = 2(Att)- d / 2 /T(D/2). 
The thermal average of any observable A is (A), the 
disorder average is A, and the average wit h th e dy- 
namical action S[u, u] is denoted by (A)s = (A). The 
Fourier transform of a function h rt of (r, t) is h quJ = 
f , e~ ig ' r+wt h r t where f. = Jdrdt, and the inver- 
sion reads h rt = J^e^-^h^, where J q = J 
f = fl 2 , The Fourier transform of the correlator 
A(u) is A K = j du e~ tK " w A(u) in general or A K = 
J a due- iKU A(u) in the periodic case. One has thus 
A(u ) = / Ke «-A K , where f K = J £ or i £ K in the 
periodic case. Note that A K is a real and even function 
of K. 



U 
U 



u 
u 



J 1 



FIG. 10. Conventions for the diagrammatics. The correla- 
tion function uu, which vanishes at T = is a full line with 
no arrow. The reponse uiu is a full oriented line. The vertex 
— ^iu r tiu rt i A(u r t — Urt' + v {t ~ O) ^ s naturally splitted in 
two half vertices corresponding to the points (r,t) and (r, t'), 
and the dashed line means that both points have the same 
position. At T = the correlation vanish. 



APPENDIX B: PERTURBATION THEORY 

We derive here the direct perturbation theory at T > 
without the use of the MSR formalism. To organize the 
perturbation series, let us multiply the non-linear part of 
the equation of motion F(r, vt + u r t) by a fictious small 
parameter a, which will be fixed to one at the end of the 
calculation. Directly on 



{Urt) 







(rjd t - cW 2 )u rt = aF(r, vt + u rt ) + f + ( r t + Kt 



(Bl) 



we can formally expand u — X) n >o oi n u^ n \ f — rjv = f = 
X)«>o a ™/ > s °l ve recursively the system (Bl), even 



at non-zero temperature, and compute the a-expansion 
of every observable. Note that we added a source h rt 
(with no constant uniform part) so as to compute the 
response function. As the force is gaussian, the expan- 
sion of disorder averaged quantities is in powers of a 2 , 
and is in fact an expansion in powers of A. We denote 
by C r - r ' t t-t' = (urtUr't/) the exact correlation and by 
r - r i_t-t' — ( sh*','*, ) ^ e exac t response functions. 
The first iterative steps are = = and 



n 





-cV 2 


4? 


— Crt + h rt 




-cV 2 




= F(r, vt + i$) 


{r]d t 


-cV 2 




= d u F(r,vt + u ( ° t > )u { r 1 t ) +.f (2 



These are sufficient to compute to first order in A the 
force, the correlation and response. 

In the absence of disorder the system moves with a 
linear characteristics / = ryv and one has the following 
correlation and response 
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c, 



2r]T 



" ' (cq 2 ) 2 + (r]uj) 2 
R - 1 

-L^quj — 9 

cq — irju) 



-cq \t\/ri 
C « = T 13 ( B2 ) 



cq* 



R 



qt 



V 



-cq t/r) 



(B3) 



related by the fluctuation-dissipation theorem (FDT) 
TR rt = —9(t)dtC r t- Note that 1Z and C do not verify 
FDT at v > 0. 

To first order in A one obtains at T = 

f~W = - I (B4) 



K<3 eg- - 
^ ~ (cq 2 ) 2 + ( V lj) 2 



V = 



1 



q '" (cq 2 ) 2 + ( V U) 2 



(B5) 

A"(vt)Ro t (l - e*"*) (B6) 



These results can be extended to any temperature T: 

/-t^ = -X>i(o; = 0) (B7) 

(B8) 

Rq, u V ((j)R-q, 

#g,u, (X> 2 (w = 0) - Da(w)) C g , w + fc.c. 



(B9) 



(2>a(w = 0)-Z>aH)i2 g , 



where we have introduced the effective vertices, non-local 
in time, smoothed by the temperature (see Figure [ll]) 

V (t) = ( A K e iKVt+< - lK)2{Co °- Cot) (BfO) 

J K 

2?x(t) = J i K A K e iKVt+ilK)2{Co °- Cot ' ) R ot (Bff) 

J hi 

V 2 (t) = { (iK) 2 A K e iKVt+< - tK)2{Co °- Cot) Rot (B12) 

J K 

— = + 



+ 



+ 



+ 




FIG. 11. Computation of the correlation (top) and re- 
sponse (bottom) functions to first order in perturbation the- 
ory. At T > 0, the tadpoles and self-contractions of the ver- 
tices contain an arbitrary number of correlations. 



We now want to compute the corrections to the pa- 
rameters c, rj, /, T, A(u) so that v, C quJ , lZ qoJ remain 
unchanged while the physical (ultra-violet) cutoff A on 
the q integrations is reduced. To first order in A and T, 
one obtains, 



dc = 



dt} = - JtR>A"(vt) 



df 

dT - 



i? >AV) 



1 



V Jt>o 
dA(u) = C>A"(u) 



tC>d t A"(vt) 



(B13) 
(B14) 

(B15) 

(B16) 
(B17) 



with d = — Aj£ and R> t , C> t are the on-shell gaussian re- 
sponse and correlation functions, i.e., with modes q lying 
only between A — dA and A. 

A completely different way for obtaining the pertur- 
bation expansion is presented in Refs. |54|,||, as a first 
attempt to include thermal fluctuations in the large- 
velocity expansion of Ref. 18. It consists in splitting the 



displacement field into a T = part and a thermal part. 
This procedure is probably only true to first order in T 
and not controlled at higher T. Instead, the method pre- 
sented here is really an expansion in disorder at any T . 

Although the calculation can in principle be pushed to 
second order, the method is too cumbersome to do it in 
practice (see however Ref. ^ at T = 0). It is easier to 
use the formalism of dynamical field theory as shown in 
Appendix ^|. 



APPENDIX C: DERIVATION OF THE FLOW AT 
FINITE VELOCITY AND FINITE 
TEMPERATURE 



Here we give the details of the renormalization proce- 
dure used for the moving system. We u se the MSR for- 
malism with action S[u, u] given by (4.2). Having shifted 
the field u rt so that its average vanishes 



(u rt ) = 0, we 
can do perturbation theory with the gaussian part 

S [u, u] — / [iiirt (r]dt - cV 2 ) u rt - rjTiu rt iu rt \ (CI) 

Jrt 

of the action. The gaussian correlation C r t and response 
R r t functions were defined in Appendix [b|. 

The interaction part of the action contains the disorder 
correlator and also the pinning force / = 0(A): 



Si[u,u] = -f 



lUrt 



(C2) 



iu r tiu rt A(u r t — u rt > + v(t — t')) 
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The effective action for slow fields u, u is given by the fol- 
lowing cumulant expansion where the averages are com- 
puted within the gaussian part Sq over the fast fields 



m > , u > 



S<[u,u] = Sq[u,u] + (Si [u + u , u + u ]) 

-^(Si[u + u>,u + u>} 2 ) c + O(Sf) (C3) 

We now turn to the computation of the first and second 
order terms. 



1. First order 



To first order . the corrections arise from the graph 
shown in Figure [ij: 




FIG. 12. First order RG corrections. The internal lines 
carry fast fields. 



They read 

(Si[u + u > , u + u > ]) = — / / iu rt 

Jrt 

{in)A K [u] (r, t, t')R£ t _ t Au r t 

A K [u](r,t,t')iUrtiUrt< 



(C4) 



rtt' k 

l 



' rtt n 

with the shorthand notation: 



A K [u](r,t,t') = A rte «("'"-", t '+"( t -*')) e ^ ( c oo-c ot _ t ,) 

The term ( |C4| ) appears to be the sum of a iuJ-[u] term 
and a iuiuQ[u] term. Let us begin to deal with the first 
type. A short time expansion of e lK ( u rt-u rt >) y[ e \ds the 
following operators 

-( / iu rt ) f i K A K e lKVt e^ K)2{c oo- c ^R^ t (C5) 

Jrt J Kt 

which is a correction to / and 

- ( / iu rt dturt) I (in) 2 A K e {l ^^~ c ^tR> (C6) 



which is a correction to r\. The elasticity operator iiV 2 « 
is not corrected and no higher gradients like raV"a are 
generated in the equation of motion. Note aka that to 
this order, no KPZ term iu(Vu) 2 is generated^. 
The iuiuQ[u] term can be rewritten as the sum of 

-- [iu rt iu rt , [ A K e (iK)2c o > "e lK( " rt ~""*' +t ' ( *- t ' )) 

2 JrW J k, 



which has the form of a disorder correlator and yields a 
correction to A (it), and an operator quasi local in time 



lU r tlU r t* X 



(C7) 



^ e (i«) CQ Q -\-iK,(u r t— u rt /+v(t — t )) 



/v, / 1 - e o*- 



which yields a co rrection to the J . iu r t%u r t term. The 
projection of (C7) on this thermal noise operator is 



lU rt lU rt | / A K 



1 — e K c ot 



To obtain the correction to the temperature T, one 

=00 ST _ S n T _ Sr, 



T r)T r) 



An integration by parts of (|Cq), 



ST 



thanks to FDT for the "pure" R and C, yields T 
To summarize, 

Sc = 

5f= { iKA K e( iK f( c °°- c &e iKVt R$ t 



6r)=- [in) 2 A K e {lK) ( c oo-c ot ) e iKVt tR$ t 

J Kt>0 

SA(u) = I A K e (lK)2c ^e lKU 
V 5T^ [ iKvtA K e^ c °°e iRVt (l-e K2 °°A 

'Kt>0 ' 



The correction Sf has the same form as the perturbative 
expression for /, with opposite sign and shell-restricted 
functions C,R. Note that Srj = —-^Sf. 

In the infinitesimal shell limit, the shell-restricted func- 
tions C > , R > which arc evaluated at r = are of order 



dl. The differential flow is thus given by (4.11) 



2. Second order 



The fast-modes average (S 2 ) c can be decomposed into 

one term with / in factor plus the rest which does not 
contain /. The former vanishes for the following rea- 
son: the contraction of the / J t iu r t with the u r i tl or 
u r 't 2 contained in the vertex operator involves a fast re- 
sponse R>._ rt . But Jl-Rrt — 0, since its modes live 
in the shell. The latter is the connected average of 
two disorder vertices. We now extract from it a cor- 
rection to the disorder, i.e., a term which has the form 
— \ J rtt , iurtiurt' 8 A(u r t — u r t> + v(t — t')). The corre- 
sponding diagrams are represented in Figure [13l 
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V 



a 



b 





FIG. 13. Diagrams needed to compute the second or- 
der corrections at T = 0, and at any T in Wilson's scheme. 
Each of the two external lines (corresponding to a u field) is 
connected to a tree of response functions (the lines) due to 
causality, and provides an orientation to these lines: we drew 
the arrows just for clarity. 

Each diagram has two external iii r t iu r >t' legs, to which 
corresponds a functional half vertex of u r t and u r 't' re- 
spectively. Calling r, t' the (positive) time arguments 
of both response functions, denoting U = u rt — u r i t > + 
v(t — t'), the diagrams have the following analytical ex- 
pressions, integrated over r, r', t, t', p, r, r'; 

a = -iu rt iu rV 5 r ^ r ,A"{U)A(U + v{r' - r))i?> R> T , 
b = -iu rt iu rV 5 r >- r - p A'(U + vt')A'(U - vr)R> T R> pT , 

C = iu rt iUr't'Sr'-rA"(U)A(v(T' - t))R> t R> t , 

d = iu rt iu r , t ,5 r ,- r - p A'{U + vt')A'(-v(t' + t))R> t R> pt , 

After another short distance expansion of b and 
d, noting that J p R^ pT R> T , = S q R% qT R> T > = 
SoA- D e~ cA ( T +T '/' 7 (iZ /?7 2 , a pr oper symmetry counting 
yields the term of order A 2 of ( 4.11 ). 

The results obtained here are consistent with the anal- 



ysis of Rcf. 56 



APPENDIX D: NEW RESULTS IN THE 
NON-DRIVEN CASE AT FINITE 
TEMPERATURE 



We give here a detailed analysis of the functional renor- 
malization group flow at T > and zero velocity. The 
temperature is an irrelevant operator and flows exponen- 
tially fast to zero. We show however that the temper- 
ature rounds the cusp in a region of size proportional 
to T around the origin and that in this boundary layer, 
the disorder correlator takes a super-universal (to low- 
est order in e) scaling form. In addition we show how 
to carry a systematic expansion at low T . As temper- 
ature decreases, the correlator of the disorder becomes 



more and more pinched, and eventually reaches its zero- 
temperature cuspy fixed point at infinity. 

We show that during the renormalization at v = 
with a flowing temperature Ti — * 0, the cusp forms only 
asymptotically (I — > oo), and A(u) has the following scal- 
ing form in the boundary layer \u\ ~ Ti/x 



A,(u)~ Ai(0)-T,/(ux/Tj) 



(Dl) 



with f(x) — y/1 + x 2 — 1 and where \ = |A*'(0 + )| mea- 
sures the cusp. 

Furthermore, we show that the following expansion in 
temperature for the solution of the FRG flow holds 

{A{u)-A{Q)-Tf = Y,T n Uu/T) (D2) 

n>2 

thus we obtained a fairly complete picture of the solution. 

1. The curvature 

The flow equation of the value at zero of the disorder 
correlator is 



d l A l (Q) = {e-20A l (0)+T l A' l '(0) 



(D3) 



Since A/ — > A*, the convergence of A/(0) towards A*(0) 
implies that T; A"(0) also converges. From the fixed point 
equation 

(e - 3C)A» + C («A*(tO)' = ~ (A*(«) - A*(0)) 2 " 

one has simply (e - 2C)A*(0) = A*'(0+) 2 , and thus, 

-TiAf(O) -> A*'(0+) 2 (D4) 



2. Scaling function in the boundary layer 

We show here that the assumption that the curvature 
at zero of A; diverges like a power of the inverse temper- 
ature implies that all the derivatives at zero also diverge 
and that there exists a well defined and particularly sim- 
ple scaling function in the boundary layer around zero. 

Precisely, for any function Tj decreasing to zero and a 
function A; (u) such that 

fyAtiu) = (e - 2C)A ; (n) + (nA^u) + TiA"(u) (D5) 

+ a;» (a ; (o)-a ; ( u ))-a;( u ) 2 

if AJ'(O) - -{ X 2 /Ti) a for some a > 
and x, then, defining the functions fi{x) = 

jr (Ai(0) - A l (xT { l a+1)/2 x~ a )), we obtain that every 
derivatives of // at x = converge to the corresponding 
derivatives of f(x) = \/l + x 2 — 1, and that / is the only 
fixed possible fixed point for /). 
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A simple way to see the convergence to the scaling 
function / is to write the flow of // 



1 



T 1 A( / (0) + -x ao T I 1 - Q (l + / I ) 



Ti difi - 2/, + 



■in 



and eliminate at large / the rhs term which is subdom- 
inant (higher order in Tj) for a > 0, since T has been 
absorbed in the variable x of /;. We have used that 
6 = 2 — e + 2£. Hence the fixed point equation for /; is 

|(1 + /) 2 " = 1 

which has the solution f(x) above since we know that 
/(0) = 0, /"(0) = 1 and / (4) (0) = -3 is easily checked. 

This is confirmed by the study of the flow equations 
for the successive derivatives a n — A^ 2n ^(0): 



da n = (e + 2(n - l)()a n 

+ Ta n+1 - 5 J2k=i 



(D6) 



2(nH 
2k 



a fcOn+l-fc 



From a trivial recurrence, the hypothesis A"(0) 
— (x 2 /T;)" implies that T™( Q+1 )~ 1 a„ converges for any 
n. Moreover the limit c n 
can be obtained from ( |D6| ) and is c 

To fix the value of a (a = 1 as strongly suggested by 



lim^ 00 T"( tt+1 )- 1 X - 2 ™ Q a„ 

_ (1.3...(2n-l)) 2 _ 
2n-l 



(D4)), we checked that the only values of (3 > 0, 7 > 
such that gi(x) = ^ r (A*(T l ^x) - Aj(lf a;)) has a mean- 
ingful fixed point are (/3, 7) = (1,1). For these values, 
the fixed point is g(x) = A*'(0+)x + Jl + (A*'(0+)x) 2 . 



3. Next order in T 

The procedure which gives us the leading behavior in 
the boundary layer controlled by temperature can be ex- 
tended analytically with arbitrary accuracy in an expan- 
sion to any order in T. We study 



diAi(u) - (e-2C)A,(u) + C«A{(«) 



a;» (a,( u )a j (o))-a;h 



with 9 = 2 — e + 2£. For numerical purposes or 
for the following analytical computation, it is useful 
to switch to the function y(u) = (A(u) - A(0) - T) 2 
which remains quadratic at the origin when T — > 0, 
since y(u) = T 2 + |TA"(0)|u 2 + 0(u 4 ) for T > and 
y(u) = A'(0+) 2 u 2 + 0(u 4 ) for T = 0. This function 
flows as 

8 iy = 2(e - 2()y + Quy 1 +Jy{y" - y"(0) - AT) 



We can replace the scale I dependence of yi(u) by a T 
dependence since T and I are linked by T; = T§e~ el . The 
function yr(u) can be expanded in 



y T (v) = Y,T n f n &) 



T 

The expansion begins at n = 2 since yr(0) = T 2 and we 
have 

/a(0) = 1 / w >a(0) = 
The equation for the /„'s reads 

T n ((2(6 - 2C) + n0) f n + (C - ejar/O = 



/£T»/n 4t~^t»- 2 (/:-/;'(o)) 



One can solve this equation order by order in T. It is 
useful to divide A and T by \ 2 an d w by \- With these 
rescaled quantities, we have simply 

^(0)=-2T ; Af(0)^2 

and thus /^'(O) = 2. If we knew the full behavior of 
y!^(0), i.e., the /^'(0)'s, we could completely solve the 
system. Here, we get 



/ 2 (X) = 1+.T 2 



Mx) =41- 



1 - — 



-(4 — (e — C))a;(asinh x — x) 



c 



+/3'(0)y 

where we wrote /s(a;) such that the three first lines are 
functions which vanish and have zero curvature at zero. 
Note that while /^(x) is universal, the last term fs(x) 
contains un unknown integration constant /^'(O) which 
presumably depends on the initial condition of the flow 
and is thus not universal. Indeed we observed a non- 
universal /3(C)) in a numerical integration of the flow of 
y T (u). 

The procedure can be carried to any order in T and 
the all the / n 's are accessible. The unknown coefficients 
of the expansion 

-2T I A{ , (0) = 2 + ^T»/« +2 (0) 

n>0 

are similarly non-universal. 

Both Subsections D 1 and p 2| thus provide a rather 



convincing and consistent picture for the solution of the 
T > 0, v = FRG equations (awaiting a mathematical 
proof). 
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APPENDIX E: ANALYTICAL SOLUTIONS AT 
FIXED TEMPERATURE 

We present here the analytical solutions of the fixed 
point equations for RF and RP at fixed T. Thanks to 
the exact expression of these fixed points, we are able 
to check the scaling form derived in Appendix |d| within 
an "adiabatic" hypothesis where the running correlator 
at I is identified with the fixed point at T = TJ. Our 
families of fixed temperature fixed points (FTFP) giv 
back the known fixed points at T — in both the RFl 
and the RPU cases. However, even if we obtain the same 
form @ for the RF T = fixed point as in Ref. §, we 
disagree with the scaling in e. 



1. Random field 

We look for a fixed point of 

d^Uu) = (e - 3C)A ; (u) + C («A,(«))' (El) 
- \ (Aj(u) - Aj(0) - T) 2 " (E2) 

with fixed T and initial random field condition J Aq > 0. 
Since d\ In J du A;(u) = e — 3£, a meaningful fixed point 
can be obtained only for £ = e/3. Fixing the RF strength 
J Ao to one, we are led to the following problem: for any 
T > 0, find the fixed temperature fixed point function 
(FTFP) A(T, u) such that 

|(uA(u))' =i(A( u )-A(0)-T) 2 " (E3) 



duA(u) = 1 



(E4) 



Integrating (E3) from to oo yields TA'(0 + ) = 0, hence 
the FTFP has a cusp for T = and no cusp for T =/= 0. 

At T = 0, integrating (|E^) from to it and dividing 
by A(u) yields u = A'{u) - A(0)A'(u)/A(u). Then, in- 
tegrating again from to u yields the T = FTFP, by 
imposing (E4) 



A(T = 0,u) 




1/3N 



where the function y{x) is implicitely defined b 



— =y-l-hxy 



Since y(Q) = 1 one has 
A(0,0) = 



1/3 



3(2V2/ Vy-i-lny) 2 



(E5) 



(E6) 



(E7) 



It is easy to compute the number J y = dxy(x) = 
2 Jq 1 dyx(y) = 2^2 ft dy^y- 1 -lny ~ 1.55. Note the 



behavior near given by y(x) = 1— |x| + ^- — +C(a; 4 ) 
thus (9^A(0,0 + ) = ^. Note also the gaussian decrease 

of correlations at infinity y(x) ~ e _1_x . 

An intriguing fact is the scaling of the T = fixed 
point with e: its n th derivative at + scales like 



^A(T = 0,M = + )-e( 1+ ")/ 3 

At T > 0, there is no cusp (<9 U A(T > 0, 0+) 
the same double integration of (|E3|) yields 



A(T, u) = A(T, 0)y(T, Wv / £ /(3A(0))) 
with y(T,x) implicitely defined by 

T 



x 

Y 



= y-l-(l 



A(0) 



\ny 



(E8) 
0) and 

(E9) 
(E10) 



The value of A(T, 0) is determined by condition ([EJ). 
Using J dxy(x) — J dyx(y), this condition reads 



24A(T, 0) 3 




y-i-(i 



A(T,0) 



)lny = 1 (Ell) 



This equation admits a unique solution A(T, 0) > for 
any T > 0. Then there exists a unique FTFP A(T^u) 
for each T > 0. Som e of them are displayed in Figure |l4|. 
Note that T = in plO| , pTT| ) gives back the T — non- 
analytic fixed point A(T = 0,w) (|E5|). Hence the set of 
FTFP has a nice T — > limit, even if there is a qualita- 
tive difference between the cuspy T — FTFP and the 
analytic T > FTFPs. 

As is obvious from their analytical expression, or from 
Figure 0, limr^o A(T>) = A(Q,u): the T — non- 
analytic fixed point is approached smoothly by the set of 
analytic fixed T fixed points. When T approaches zero, 
the curvature of the FTFPs at the origin goes to — oo like 



lim -TA"(T,0) = -A(0,0) 

T^O 3 



(E12) 



with A(0, 0) given by 

We also checked that the A(T, u) converge when T — > 
to the zero t emp erature fixed point with the predicted 
scaling form (Dl) 



A(r,o)-A(r,rf; 

T 



(E13) 



where \ = |<9„A(0,0+)| is given by the T = FTFP 
equation \ 2 = f A(0,0). 

Some of the RF fixed T fixed points are shown on the 
right bottom quarter of Figure [3, including the cuspy 
(highest) T — fixed point. Absorbing e in T and A, we 
chose to plot the non-trivial solution to the most reduced 
problem 



(uA(u))' = i (A(u) - A(0) - T) 2 " 
duA(u) = 1 



(E14) 
(E15) 
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To restore e and J A, one simply has to note that the 

"dimensions" are T = A = e 1 / 3 (/ A) 2 ^ 3 and u = 

e~ 1/3 (f A) 1/3 . The left bottom of Figure || shows 
— TA"(T, 0) as a function of T. This combination has 
a finite limit (~ 0.17) when T -> 0. 



2. Random periodic 



In the Random Periodic case, the conservation of the 
period a of A requires £ = 0. After a suitable rescaling, 
u -> u/a, A -> A/(ea 2 ) and T -> T/(ea 2 ), the fixed 
point equation reads for the 1-periodic function A(u) 



A( U ) = 1(A( U )-A(0)-T) 2 " 



(E16) 



and is easily solved by quadrature, by analogy with a 
particle's position X(u) = (A(u) — A(0) — T) 2 at time u 
in a potential F(J*Q = 4A 3 / 2 /3 - 2(A(0) + T)X verify- 
ing A" (it) = — V'(X(u)). The quadrature leads to the 
reciprocical function u(X), parametrized by A(0) and T, 
as a sum of two elliptic functions. Then, imposing the 
solution A(u) be 1-periodic fixes A(0) as a function of T. 
The result is: 

• for T > (2tt)" 2 the only solution is A(u) = 

• for < T < (2tt)~ 2 another solution arises, which 
resembles a cosinus function of linearly vanishing 
amplitude when T — > (27r)~ 2 . This non-trivial so- 
lution has no cusp but becomes pinched as T de- 
creases (growing curvature |A"(0)| and higher har- 
monics). As can be seen on the analytical expres- 

sion (not given here), — TA"(0) — > 1/36. In par- 
ticular, it remains finite when the temperature van- 
ishes. 

• eventually for T — > 0, the non-trivial solution uni- 
formly tends to the zero temperature fixed point 



A(u) = -(- - u(l - u)) 
■ ' 6 V 6 v " 



(E17) 



The temperature (27r) -2 in our units is exactly tho,crit- 
ical temperature T g of the random-field XY-modelEll and 
the fixed points near T~ reproduce the line of fixed points 
of this problem (since we worked to second order, it is 
only an approximation). Indeed in D = 2, the naive di- 
mension of the temperature is zero and our FTFP has a 
direct physical meaning. Note that another random gra- 
dient term becomes relevant-,in D = 2 but does not feed 
back on the flow of A(u)E3G3. 




-TA"(0) 
-0.17 





1 



1 



FIG. 14. Fixed points at fixed temperature A(T,u). Bot- 
tom: RF case, the solution A(T, u) to the reduced equation 
( B14 , ^15| ) exists for any T. Right bottom: plot of the RF 
FTFPs for T e {0,0.1,0.2,0.5,1}, these temperature are 
located on the abscissas of the left bottom plot. On the 
left, plot of -TdlA(T,u = 0) versus T. Top: RP case, 
the 1-periodic non-trivial solution A(T,u) to ( |E16| ) exists for 
< T < (2tt)- 2 . Right bottom: plot of the RP FTFPs for 
T e {0, 0.005, 0.01, 0.015, 0.02}, these temperature are located 
on the abscissas of the left top plot. The FTFPs are analytic 
for T > but tend to the cuspy A(T = 0, it) as T ~* 0. The 
curve on the left shows —Td^A(T, u = 0) as a function of T. 
It is not a straight line. 

We can now use the exact FTFPs to check that an 
adiabatic hypothesis is consistent with the scaling form 
(Dl). Indeed, one can numerically check that the cor- 
relator with a flowing temperature has the FTFPs have 
the scaling (Dl) as T — > 0. To conclude about the prob- 
lem with a flowing temperature T; — > 0, it appears from 
these observations that no cusp occurs at finite scale for 
To > 0. The cusp forms only asymptotically (I — > oo), 
with 

lim -TiA{'(0)=A(0,0) 

/ — >OG 



given by (E7) and it obeys a scaling form in the boundary 
layer |u| < Ti/x 



AtW-AiW-Tifiux/Ti), f(x) = V^ + x 2 -l 
We note that the precise form of the flow of the tem- 
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perature (i.e. the value of 8) only affects sub domi nant 
behavior (i.e. the function /a (a:) in Subsection (D 3)). 



APPENDIX F: MULTI-DIMENSIONAL CASE 

We give here a possible extension of the FRG to a 
multi-dimensional displacement field. This study gener- 
alizes the approach of Ref. |5^ by including the effect of 
v > and T > in the flow. For periodic structures, a 
similar study of the multi-dimensional displacement field 
was shown in Ref. [56| to yield novel effects. 

In a D + N dimensional space, we distinguish between 
the internal or longitudinal space of dimension D, to 
which r belongs, and the transverse space of dimension 
N, to which u belongs. The elastic energy of an inter- 
face without overhangs defined by a height function u r is 
quadratic in Vu of the form ~ J r c^" ' {d tl u t r ){d v ui r ). 

The disorder: the random bond (RB) case corre- 
sponds to a random potential V(r,u), with correlations 
(V{r, u) - U(r',it')) 2 = -25° (r - r')R(u - u'). Func- 
tion R(it) is even, vanishes at u — and goes to a neg- 
ative constant for \u\ 3> r/. The random field (RF) 
case corresponds to a force F l (r, u) with correlations 
F^r, u)Fi{r',u') = 8 D {r - r')A ij (u - u'), where the 
A lJ (it) vanish for |u| ^> r/. A RB gives rise to a 
RF via F l = —d l V and the correlators are related by 
A ,j (m) = — 9 y R(u). Note that this type of correlator 
deriving from a RB has J d^uA'^ti) = 0. Finally, 
the random periodic case (RP) occurs when it is de- 
fined up to a discrete set of translations forming a lat- 
tice of points P, e.g. when it is a phase, defined up to 
2ir shifts. In this case, the disorder is periodic and one 
has A ij (u) = A lj (u + P) for any P of the lattice (or 
R(u) = R(u + P)). 

The overdamped dynamics is given by 



rfjd t u 



o _ 

rt 



ii t + F l (r,u rt ) + Crt + f l + K 



where r\ is the friction tensor and C a Langevin noise, 
with correlations (QtCt') = 2(r)T)^5(r - r')S(t - t'). 
The tensor T stands for the temperature(s) of this out 
of equilibrium system. We added a driving force f l per- 
pendicular to the interface and a source field h l rt , as an 
external excitation. 

Without assuming any symmetry, let Cl\ and R l r 3 t be 
the gaussian correlation and response functions. We ob- 
tain by the same procedure as for the N = 1 case the 
following first order corrections due to disorder 

Sc lj = 

J Kt 

Sr) ij = _ j e ^(Coo-c ot )A K +i K .vt A ^ iK i tR ik iK j 

gire.tit/gire. (Coo-Cot). in _ e iK.C 00 .iK\£ij 



Using (A™' 3 )* = Af£, All = we write the on- 

shell corrections as (with the matrix product A.B = 



K» = 




*/=- 


/ in. 




J K 


5i] = — 


/ in. 




I K 


S(2 V .T) = - 


/ in. 




J K 


SA K = in 


Cg .i 



D> A JlK.Vt 



l>.A. K e iK - vt i 



1> tv-A J-K.vt 



The second order correction to A reads 



6A lJ (u) 



qrr' 



r> > rn k. jD>7n'l 
H qr K qr' 



5(T)T) 



(A kl (v(T' - t)) ~ A kl (u + v(t' - t))) d m d m 'A ij (u) 

-d m A il (u + vT)d m ' A k] {u - vr') 
-d m A u (u + vr)d m ' A ]k {v{T' + t)) 
+d m A a {v(T' + T ))d m ' A kj {u - vt') 



Note that each of the first three terms are symmetric 
under i <-» j, u <-> — u and that the fourth is exchanged 
with the fifth under this symmetry. Then A remains a 
correlator. 

Of course this second order correction to A gives back 
the expression already computed for a D + 1 interface if 
N = 1. At zero velocity, one gets the seaond derivative 
of the flow equation of Balents and Fisherfl If we assume 
that A (it) depends only on the component of u parallel 
to the velocity and send v to zero, then our .expression 
reduces to the equations of Ertas and KardarEj. 

To simplify the analysis, let us rely on the assumed 
symmetries of the system. If we suppose that the initial 
problem is rotationnally invariant, i.e. has O(N) sym- 
metry, then the elasticity tensor c, the friction tensor 77 
and the temperature tensor T are only scalars and the 
force-force correlator A is covariant, i.e., for any 1Z such 
that K^.K = 1, TZ^.A(u).n = A(K.u). 

During the flow, we expect from physical grounds that 
the running terms of the action will conserve their sym- 
metries but the velocity v which is fixed once for all se- 
lects a particular direction in transverse space. The in- 
teresting symmetries are given by the little group of the 
velocity, i.e., the transformations TZ such that TV .TZ = 1 
and IZ.v = v. Then one may decompose the tensors on 
a basis involving v (one has only two frictions, tempera- 
tures, response and correlation functions and fivet2l Aj's, 
functions of (it 2 , v 2 , u.v)). 

Unfortunately, the full problem can not be easily de- 
coupled, even with the simplifications pointed out above. 
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No closed equation e.g. for the correlator restricted to 
displacements aligned with the velocity A(u \\ v), has 
been found, and the problem even at zero temperature 
seems involved. The simplification used in Ref. ^ con- 
sists in assuming that A does not depend on the trans- 
verse coordinates. This assumption reduces the problem 
to the N = 1 case, and it would be interesting to solve 
at finite T the behavior of transverse coordinates along 
the lines of our analysis. 



= eA"(0)-3A"(0) 2 

-9A 2 A"(0)A iv (0) + 0(A 4 ) 

(note that A"(0) < whereas A iv (0) > 0). The flow of 
the friction is similarly slowed down: 



d In r\ 



[ e-°(2-s 

Js>0 



AjXs) - A(0) 
A 2 



A"(O)-3A 2 A IV (O) + 0(A 4 ) 



APPENDIX G: THE FLOW OF THE DISORDER 
CORRELATOR AT SMALL VELOCITY 

The effect of a small velocity on the FRG flow is mainly 
restricted to the boundary layer of width pi about the 
origin. Analytically, it is rather difficult to give an esti- 
mate of pi or to decide how A/ (u) precisely behaves in the 
boundary layer |u| ~ /?;. It is however possible to simplify 
the formidable second order correction to the disorder 



correlator, displayed in (4.11), and to obtain analytically 
several results, giving s ome hints about this behavior. 

The C(A 2 ) term in (4.11) is written under a form in- 
volving two integrations over s, s', reflecting the presence 
of two response functions integrated over time. After 
some integrations by part, the 0(A 2 ) term becomesE] 

A»/e-* (A(Xs) - A( " + AS) + A( "- AS) ) (Gl) 
c _ s A(u + \s)-A(u) j- e _ s A(u-\ s )-A(u) 



s>0 ^ Js>0 ^ 

c - s A'(\c) A{u + Xs) + A(U " Xs) ~ 2A(u) 

s>0 



A 



Integrated over u, this correction becomes 

A(u + As) - A(u- As) 



duA'(u) / e~ s (2-s)- 



s>0 



A 



For any non-crazy function A, this expression is positive. 
Assuming that A has no cusp, it can be safely expanded 
and we can check that it is of order A 2 : 



d / A = (e - 3C) / A + 2A 



A' 



G(A 4 ) 



Thus at v 0, the integral of A grows during the flow, 
whereas it was conserved in the statics. 



Using (Gl), one can also compare the flow of A; at 
small velocity to the cuspy v = flow. In particular, one 
observes that the effect of the velocity is to reduce the 
blow-up of the curvature A"(0) 



3A"(0) = eA"(0) - A"(0) / e" 

Js>Q 

'A'(As) V 



, A(A a )-A(0) 
A 2 



e 

s>0 



A 



The term -A'(0+) 2 in the v = flow of A(0) in (|4.14| ) 
at v = is replaced at v > by 



9A(0) = (e - 2C)A(0) 




.- s A(As) 



s>0 ^ I Js>0 \ ^ 



A(As) 



(G2) 

2 



which has the right sign (using Cauchy inequality) to 
slow down the exponential growth of A;(0). 

Obtaining a numerical integration of the flow is a 
highly non-trivial questo, since all the interesting prop- 
erties occur close to the origin yielding unaccurate results 
in real space. In Fourier space, the number of harmonics 
to be retained is huge if one wants to focus on the quasi- 
cuspy behavior (A* ~ k~ 2 at the cuspy fixed point). 
However, we obtained, at least at the beginning of the 
flow (up to l c ) with small initial velocity, the curve shown 
in Figure [| The initial condition was a RB disorder (full 
line). It is obvious on the snapshot (dotted line) close to 
L that the flow transformed the RB into a RF. 



APPENDIX H: BEFORE THE LARKIN LENGTH 

We show here that at the scale l c , A; (it) is very close 
to the static zero-temperature fixed point A*(u). This 
can be checked numerically, even in the presence of a 
small temperature and small velocity. Analytically, one 
cannot obtain an exact integration of the flow, but we 
can compare A/ (u) to the known A* (u) by the following 
arguments. 

Let us take e.g. the RF case for which £ = e/3. For 



weak disorder one obtains from the integration of (4.15) 



3|Aff(0)| 



1/3 



1/3 



(HI) 



where we have used |Aq(0)| ~ A o (0) /rj ~ J A/r 3 f . 
We prove in Appendix [El] that \ = |A*'(0+)| verifies 

X ~ e 2 / 3 (/ A) and that A*(0) w e 1 / 3 (j ' A^^ \ 



Thus the range 



A*(0)/|A'*(0+)| of A* verifies 
r* » rye-«" w X /e (H2) 
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To determine the range of A; c (u) we use the fact that 
at the beginning of the flow on e can neglect the non- 
linear term in the flow equation (4.11). We are left with 
A/(u) = e^ e ~ 2 ^ l Ao(ue^ 1 ) and thus the range of A; c (u) is 
simply 



r fQc) - r f e 



-Go 



(H3) 



A comparison of ( H2) and (|H3|) shows that the two ranges 
are similar. Furthermore, in the RF case J A is conserved 
by the flow at v = 0, and thus the similarity of ranges 
shows that the shape of A; c is close to the shape of A* 
(same integral, same range). 

Similarly in the RP case it is also true that A/ c resem- 
bles A*, but now x = as can be seen on the fixed 

point A*(m) = | ( ^ — u(a — u) 
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